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EXECUTIVE  SUMMARY 


This  report  documents  work  performed  in  2015  and  2016  under  Space  and  Naval  Warfare  System  Cen¬ 
ter’s  (SSC  Pacific)  Naval  Innovative  Science  and  Engineering  (NISE)  Program  supporting  the  Office  of 
Naval  Research  Expeditionary  Warfare  and  Combating  Terrorism  Department  (ONR  30)  in  its  Tactical 
Networking  and  Small  Unit  Technologies  Thrust  Area. 

This  project  quantities  performance  of  wireless  power  control  algorithms  optimizing  a  mobile  network 
in  a  ground-based,  small-unit  operation  sweeping  through  urban  environments.  The  wireless  network  con¬ 
sists  of  20  nodes  or  less:  most  nodes  arc  battery  limited,  most  operate  at  the  street  level,  most  arc  human 
mobile,  a  few  nodes  operate  on  rooftops,  one  or  two  nodes  may  be  a  drone  supporting  relay  and  video,  ur¬ 
ban  operations  guarantee  the  network  will  fragment.  These  wireless  networks  must  adapt  to  the  changing 
RF  environment  to  meet  the  following  conflicting  objectives: 

•  Maximizing  throughput 

•  Minimizing  power 

•  Guarantee  quality  of  service  (QoS). 

Performance  of  the  these  mobile  wireless  is  limited  by  their  self  interference:  increasing  power  in  one 
link  boosts  throughput  in  that  link  while  simultaneously  raising  interference  and  decreasing  throughput  in 
the  other  links.  Consequently,  Wireless  Power  Control  (WPC)  is  a  fundamental  control  mode  enabling  mo¬ 
bile  networks  to  trade  power,  throughput,  and  QoS  [37],  [13],  [49],  [72].  There  arc  other  control  modes: 
Dynamic  Spectral  Access  (DSA)  manages  the  power-rate  trade-offs  by  adaptively  selecting  a  best  channel 
[54],  and  pattern  steering  selects  antenna  patterns  maximizing  power-rate  performance  [16].  Joint  opti¬ 
mization  over  combinations  of  these  modes  is  possible:  power  and  channel,  or  power  and  patterns.  How¬ 
ever,  all  these  adaptive  controllers  must  operate  over  a  distributed  network — centralized  control  is  imprac¬ 
tical  for  small-unit  operations.  The  proposed  controllers  must  also  be  compatible  with  legacy  radios.  Con¬ 
sequently,  power-rate  controllers  are  limited  to  software  algorithmsf  and  control  modes  must  be  enabled 
by  existing  radios  or  by  bolt-on  hardware.  Therefore,  distributed  WPC  algorithms  arc  the  starting  point  of 
this  project.  Assessment  of  several  distributed  WPC  algorithms  operating  in  a  credible  tactical  scenario  is 
undertaken  with  two  objectives: 

WPC-1  Best  performance  bounds 
WPC-2  Ranking  WPC  algorithms. 

Best  performance  is  the  power-rate  performance  delivered  by  the  perfect  and  centralized  WPC  con¬ 
troller.  Our  findings  are  that  the  distributed  WPC  algorithms  use  roughly  the  same  power  as  the  ideal  WPC 
and  still  deliver  at  half  the  throughput — assuming  system  delay  is  small — increasing  system  delay  causes 
these  distributed  WPC  algorithms  to  become  unstable.  Indeed,  “all  adaptive  algorithms  can  be  broken” 
if  the  input  measurements  exceed  the  tracking  rate  of  the  algorithm.  Our  key  finding  is  that  a  time-delay 
compensation  scheme  stabilizes  these  WPC  algorithms.  We  conjecture  this  time-delay  compensation  is 
necessary  for  controllers  that  jointly  optimize  power  with  other  control  modes  and  is  the  basis  for  the  next 
phase  of  research. 
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1.  WHY  WIRELESS  POWER  CONTROL? 


Figure  1  illustrates  throughput  and  interference  conflicts  in  a  wireless  network.  Link  1  boosts  through¬ 
put  by  increasing  broadcast  power — thereby  boosting  interference  on  Link  2.  Link  2  increases  its  broad¬ 
cast  power  to  overcome  this  interference — thereby  boosting  interference  on  Link  1 .  The  network’s  self¬ 
jamming  and  uncontrolled  feedback  forces  the  nodes  into  a  power  race  that  ends  with  all  nodes  blasting 
at  full  power.  More  formally  [25]:  “Wireless  networks  are  fundamentally  limited  by  the  intensity  of  the 
received  signals  and  by  their  inherent  interference.” 


Figure  1 .  Transmission  and  interference  in  a  two-link  wireless  network. 

Wireless  power  control  seeks  to  manage  the  conflict  between  throughput  and  interference  while  deliv¬ 
ering  optimal  trade-offs  between  standard  network  objectives  [41]: 

•  Maximize  throughput 

•  Minimize  power 

•  Guarantee  quality  of  service  (QoS). 

The  payoff  from  wireless  power  control  is  substantial.  Douros  &  Polyzos  observe  the  following  [18]: 

•  "...  power  control  nearly  doubles  the  capacity  of  a  CDMA  system  [26]” 

•  “For  Mobile  Ad  Hoc  NETworks  (MANETS),  the  adoption  of  power  control  leads  to  an  over  50% 
improvement  to  the  energy  expended  compared  to  the  standard  IEEE  802.1 1  [60]” 

•  “[62]  estimates  the  signal-to-interference  plus  noise  ratio  (SINR)  gains  from  the  adaption  of  power 
control  exceed  10  dB” 

This  report  seeks  side-by-side  comparisons  of  distributed  wireless  power  control  algorithms  operating 
in  tactical  scenarios  of  interest  to  the  Office  of  Naval  Research  (ONR).  The  scenarios  are  the  mobile,  ad 
hoc  networks  supporting  small-unit  ground  forces  operating  in  urban  and  rural  environment  (Section  1.2). 
These  comparisons  are  the  trade-offs  between  power  and  throughput  under  QoS  constraints.  The  payoffs 
for  ONR  are  an  “apples-to-apples”  ranking  of  these  algorithms  accompanied  by  best  possible  trade-offs 
(Section  1.3). 
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Best  possible  trade-offs  between  power  and  throughput  are  computable  under  the  assumptions  that  a 
centralized  controller  exists  with  the  following  capabilities: 

•  Instantaneous  and  perfect  measurements  of  the  noise  and  interference  at  each  receiver 

•  Instantaneous  and  perfect  control  of  each  transmitter 

•  Homogeneous  network  traffic. 

The  best  bounds  benchmark  the  effectiveness  of  wireless  power  control  algorithms  to  compensate  the 
time-varying  mobile  channels.  The  ad-hoc  nature  of  mobile  wireless  networks  means  that  mobile  nodes 
may  drop  out,  the  wireless  channels  are  subject  to  fast  and  slow  fading,  the  signal-to-interference  and  noise 
ratios  (SINR)  can  only  be  estimated  at  each  node,  and  the  network  traffic  is  not  homogeneous.  Thus,  wire¬ 
less  power  control  algorithms  must  maximize  throughput  and  minimize  energy  while  operating  from  noisy, 
time-delayed  measurements,  noisy,  time-delayed  power  control  messages,  and  time-varying  network  traffic 
demands. 


1.1  MESSAGING  AND  LATENCY 


Wireless  power  control  requires  messaging  for  power  control.  Figure  2  compares  centralized  against 
distributed  messaging  for  a  two-link  wireless  network.  Centralized  control  requires  that  receivers  broad¬ 
cast  their  estimated  quality  of  service  to  a  central  node.  The  central  node  computes  an  optimal  power-rate 
trade-off  for  the  entire  network  and  broadcasts  the  new  power  commands  to  the  transmitters.  Centralized 
control  is  not  practical  for  mobile  wireless  networks  [13]:  “the  implementation  of  a  centralized  solution 
is  not  plausible  in  practice  due  to  the  signaling  overhead  and  high  computational  complexity.”  Distributed 
control  is  the  focus  of  this  report  and  only  requires  the  feedback  from  the  receiver  to  its  transmitter. 
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Figure  2.  Wireless  power  controllers. 


Upper  and  lower  bounds  on  distributed  power  control  are  determined  by  maximal  centralized  and  local 
control.  Maximal  centralized  control  is  equivalent  to  all-knowing  network  “genie”  having  access  to  perfect 
measurements  and  instantaneously  selects  best  possible  power-rate  trade-offs.  Local  control  typically  as¬ 
signs  transmit  power  before  the  deployment  and  thereby  provides  a  lower  bound  on  power-rate  trade-offs. 
Distributed  algorithms  must  outperform  local  control  to  justify  the  additional  overhead.  Those  distributed 
algorithms  that  come  close  the  centralized  control  arc  the  candidates  for  development. 
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1 .2  TACTICAL  SCENARIOS 


“Tactical”  is  in  the  title  of  this  report.  Tactical  mobile  wireless  networks  serving  small-unit  ground 
forces  operate  with  different  requirements  than  commercial  systems: 

•  Radios  are  single  transceivers 

•  Hai'd ware  upgrades  are  unlikely 

•  Frequencies  differ  from  commercial  systems 

•  Avoiding  detection  may  force  extereme  noise  levels 

•  Battery  power  is  always  a  problem 

•  Jamming  will  be  intelligent  and  hostile 

•  Nodes  routinely  drop  out 

•  Nessaging  priorities  vary 

•  Throughput  requirements  are  increasing 

•  Centralized  control  is  impractical. 

These  limitations  point  to  software-only  solutions  operating  on  the  existing  radio  hardware.  Assessing  per¬ 
formance  of  these  algorithms  in  tactical  RF  scenarios  is  a  basic  task  to  quantify  their  performance.  Power 
and  throughput  are  basic  tactical  network  performance  metrics.  Well-designed  wireless  control  algorithms 
lead  to  optimal  power-rate  trade-offs  [37].  Best  possible  bounds  on  the  power-rate  trade-offs  in  tactical  RF 
scenarios  provide  absolute  benchmarks  to  rank  system  performance  and  guide  development  in  these  chal¬ 
lenging  RF  scenarios. 

1.3  POWER-RATE  TRADE-OFFS 

The  geometry  of  the  power-rate  image  shows  best  possible  trade-offs  between  throughput  and  power 
that  any  wireless  power  control  algorithm  can  deliver  in  a  given  tactical  RF  scenario.  These  optimal 
power-rate  trade-offs  are  computed  for  a  two-link  tactical  scenario  to  benchmark  distributed  wireless 
power  control  algorithms  operating  with  delayed  and  noisy  measurements. 

Figure  3  shows  the  power-rate  image  for  a  two-link  wireless  network  operating  over  a  noisy  narrow- 
band  channel  at  300  MHz  subject  to  QoS  constraints.  This  image  shows  all  possible  power  and  rates  that 
an  “all-knowing”  network  genie  delivers  while  meeting  QoS.  Here,  “all-knowing”  means  the  genie  has 
perfect  channel  measurements  and  implements  power  control  without  delays.  The  genie  sets  the  powers  p\ 
and  p2  for  the  two-link  wireless  network  to  meet  the  QoS.  The  QoS  constraint  requires  the  throughputs  on 
both  links  to  exceed  a  minimum: 

Rb,i,Rb,2>  f-58  [Mbps], 

Equivalently,  the  SINR  both  links  must  exceed  3  dB.  The  horizontal  axis  is  the  total  network  power:  py,  = 
Pi  +  P2-  The  vertical  axis  is  total  network  throughput:  Ry  =  R^  i  +  Rt,:2-  Section  7  details  the  explicit 
parameterization  of  this  image. 
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Figure  3.  Power-rate  image  for  a  two-link  wireless  network  meeting  QoS. 


Vertical  and  horizontal  “slices”  of  the  power-rate  image  correspond  to  power-allocation  modes.  If  the 
total  network  power  is  fixed,  the  network  genie  seeks  to  maximize  total  network  throughput.  Equivalently, 
the  optimal  power  allocation  is  the  top  of  a  vertical  “slice”  of  the  power-rate  image.  If  the  total  network 
throughput  is  fixed,  the  network  genie  seeks  to  minimize  the  total  network  power.  Equivalently,  the  opti¬ 
mal  power  allocation  is  the  left  end-point  of  horizontal  “slice”  of  the  power-rate  image.  This  geometric  ap¬ 
proach  is  not  in  the  wireless  literature — to  the  best  of  our  knowledge  [75],  [79].  Computation  of  this  Pareto 
front  supplies  the  best  bound  on  any  wireless  power-control  algorithm. 


1.4  EVOLUTION  OF  WIRELESS  POWER  CONTROL  FOR  MOBILE  NETWORKS 


The  literature  on  wireless  power  control  is  extensive:  surveys  from  the  last  two  decades  include  [32], 
[53],  [66],  [41],  [1 1].  This  section  discusses  wireless  power  control  algorithms  applicable  to  the  mobile,  ad 
hoc  networks  used  by  small-unit  ground  forces  operating  in  urban  and  rural  environments.  The  emphasis 
is  on  distributed  algorithms  operating  in  the  physical  layer — there  is  no  centralized  controller,  joint  power 
control  and  scheduling  [54],  joint  power  control  and  channel  selection  [54],  cross-layer  optimization  [42], 
[6],  cognitive  radios  [38],  or  cognitive  networking  [80]. 

The  seminal  wireless  power  control  algorithm  was  authored  in  1993  by  Gerad  J.  Foschini  and  Zoran 
Miljanic  [22]:  A  Simple  Distributed  Autonomous  Power  Control  Algorithm  and  its  Convergence.  The  Fos¬ 
chini  &  Miljanic  (FM)  algorithm  updates  the  power  pi(t )  on  each  link  (I  =  1,  . . . ,  L)  so  that  the  measured 
SINR  7 tracks  the  external  SINR  jx,i  set  by  networking  requirements: 


stepsize 

Pl(tn+l)  =  Pl(tn )  +  (1  -  P)  {lx, l  ~  lM,l(tn ))  • 
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The  overline  is  the  standard  dB  notation 


P(t)  =  10  x  log  10(p(t)). 

The  FM  algorithm  seeks  to  deliver  a  target  QoS  during  communication  along  each  link.  Descendants  of 
this  algorithm  arc  listed  in  the  first  line  of  Table  1 .  One  highlight  is  a  multiobjective  version  of  the  FM 
algorithm  that  trades  power  against  QoS  [20].  All  wireless  power  control  algorithms  suffer  from  the  fol¬ 
lowing: 

•  Noisy  measurements  in  both  the  SINR  and  in  the  control  systems 

•  Feasibility — there  may  be  no  power  allocation  delivering  the  requested  QoS 

•  Time  delays  in  both  the  measurements  and  the  control  systems. 


Noisy  Measurements'.  Idealized  systems  provide  instantaneous  measurements  of  the  SINR.  Practical  sys¬ 
tems  typically  measure  the  only  received  signal  power — the  sum  of  the  signal  power,  the  interference  and 
the  noise — to  make  indirect  estimates  the  SINR.  Shoarinejad,  Speyerm,  and  Pottie  [68]  produced  an  excel¬ 
lent  paper  showing  how  Kalman  filtering  delivers  improved  predictions  of  SINR.  This  report  assumes  ideal 
SINR  measurements. 

Feasibility :  One  solution  to  the  feasibility  problem  are  the  Opportunistic  Power  Control  algorithms  that 
adapt  the  QoS  to  a  range  rather  that  a  fixed  number  [12]:  When  the  channel  does  support  a  QoS,  the  QoS 
can  be  increased;  when  the  channel  does  not  support  the  QoS,  the  QoS  is  decreased.  Such  Kalman  filtering 
algorithms  update  the  power  pi(t)  on  each  link  so  that  the  improved  SINR  measurement  7 M,i(t)  tracks  the 
time- varying  “target”  SINR  77  ;  (t): 

Pl(tn+ 1)  =  Pl(tn )  +  (1  -  P)  (7 T,l(tn)  ~  7 M,l(*n))  • 

Selected  Kalman  filtering  schemes  are  listed  in  the  second  line  of  Table  1 .  Limitations  of  the  Kalman  filter 
are  observed  in  [84] : 

•  “. . .  power  control  schemes  using  the  Kalman  filter  are  effective  only  if  the  . . .  channel  variation  and 
measurement  noise  are  Gaussian  distributed.” 

•  “Traditional  power  control  schemes  using  power  iterations  . . .  are  not  appropriate  to  compensate 
Rayleigh  channel  fading  because  all  the  schemes  require  a  long  time  to  converge.” 

The  //-.  H°°,  and  H2 /H°°  algorithms  drop  the  Gaussian  assumption  by  assuming  channel  variations  and 
measurement  noises  are  bounded.  Selected  H°°  algorithms  are  listed  in  the  third  line  of  Table  1 . 

Time  Delay:  All  adaptive  algorithms  will  fail  if  the  updates  arrive  too  slowly  to  compensate  for  the 
changes  in  the  system.  Informally,  the  measurements  are  “stale”  or  that  the  time  delays  are  not  spanned  by 
the  system  models.  The  significant  time  delay  in  wireless  power  control  are  caused  by  the  signal  process¬ 
ing  and  messaging — not  in  the  over-the-air  travel  time.  Figure  31  details  that  the  receiver  must  measure  the 
power,  estimate  the  SINR,  undertake  a  control  computation,  message  the  power  control  command  back  to 
the  transmitter,  the  transmitter  must  receive  the  message,  extract  the  power  command,  and  reset  the  power. 
This  process  consumes  several  milliseconds.  Such  delays  will  cause  the  wireless  power  control  algorithms 
to  fail  depending  on  the  channel  correlation  and  the  update  times.  Campos-Delgado  &  Luna-Rivera  ob¬ 
serve  [7]: 
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“Therefore,  the  main  challenge  in  distributed  power  control  is  to  compensate  the  time-varying 
profiles  of  the  channel  gains  and  interactions  among  all  active  users  to  achieve  the  objective 
SINR's  and  to  counteract  the  effect  of  time  delays  induced  by  the  QoS  measurements.” 


One  approach  to  accommodate  time  delays  in  the  Kalman,  H 2,  and  the  H°°  filters  is  to  expand  the  state 
vector  [48],  [45],  [48].  Another  approach,  similar  to  the  PID  controller  [2],  models  the  time  delay  itself 
without  increasing  the  dimension  of  the  state  vector — the  paper  by  Subramanian  and  Sayed  is  a  masterful 
exposition  of  this  approach  [73]. 

Gunnarsson,  Guustafsson,  and  Blom  produced  a  superb  analysis  of  time-delay  (TD)  effects  and  offer 
an  effective  time-delay  compensation  of  the  FM  algorithm  that  takes  the  form  [29] : 

Pl(tn+ 1)  =  Pl(tn)  +  (1  -  P)  (lx, l  ~  7 M,l(tn))  , 

where  7 M,i(tn)  is  an  estimate  of  the  delayed  SINR.  Table  1  lists  this  FM-TD  algorithm  in  the  fourth  line 
along  with  other  references  specializing  Smith  control,  H°°  control,  and  Kalman  filtering  to  time-delay 
compensation. 


Table  1.  Selected  distributed  wireless  power  control  algorithms. 


1 990-1 999 

2000-2010 

2011-2015 

Algorithm 

[22],  [82],  [23],  [83] 

[20],  [1],  [12] 

[36],  [68],  [70],  [72] 
[48],  [84] 

[29],  [64], [72],  [48],  [2] 

[81] 

[40] 

[13],  [49],  [34] 
[31] 

FM 

Kalman 

h'2/h°° 

Delay  compensation 

1.5  PITFALLS  OF  ADAPTIVE  ALGORITHMS 

If  multiple  transmitters  are  updating  broadcast  power  based  on  noisy  and  delayed  signals  adapting 
to  time- varying  channels — where  the  time  variation  is  determined  by  the  speed  of  the  mobile  network — 
network-level  questions  arise  regarding  the  algorithm  performance  with  respect  to 

•  Stability 

•  Donvergence 

•  Latency  effects 

•  Measurement  errors  and  control  errors. 

The  point  is  that  every  adaptive  algorithm  can  be  broken — either  the  time- variations  of  the  channel  exceed 
the  scope  of  the  control  algorithm  or  the  measurements  arc  too  corrupted  by  noise. 

This  report  assumes  perfect  measurements  and  controls  to  focus  algorithm  performance  when  adapting 
to  the  urban  channels  and  system  delays  exceeding  the  update  rate.  The  selected  algorithms  show  rela¬ 
tively  similar  performance  at  moderate  system  delays  but  fail  upon  exceeding  a  delay  threshold.  The  pay¬ 
off  of  this  report  is  the  introduction  of  a  time-delay  compensation  scheme  that  extends  the  FM  algorithm  to 
handle  larger  delays. 
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1.6  ORGANIZATION 


Section  2  lays  out  the  signal,  interference,  and  noise  models  for  a  narrowband  wireless  network.  The 
upshot  is  that  a  wireless  power  control  is  governed  by  the  channel  gain  matrices,  the  noise  level,  and  the 
system  delays. 

Section  3  describes  the  RF  propagation  in  a  three-dimensional  urban  environment.  The  gain  matri¬ 
ces  rule  the  mobile  wireless  network  performance.  The  simulated  channels  exhibit  slow  fading  caused  by 
building-size  shadowing  and fast  fading  caused  by  multipath  interference  [68]  at  300  MHz,  1.9  GHz,  and 
2.4  GHz. 

Section  5  describes  a  wireless  power  control  system  and  its  consequent  model  of  network  performance. 
This  model  of  a  wireless  power  control  system  is  the  common  operational  framework  for  the  forthcoming 
distributed  wireless  power  control  algorithms.  This  common  operational  framework  permits  comparisons 
without  getting  entangled  in  system-specific  implementations. 

Section  6  reviews  selected  results  in  matrix  inequalities.  Section  7  uses  these  matrix  inequalities  to 
compute  power  allocations  for  a  static  channel.  Equivalently,  this  section  describes  the  power-rate  image 
shown  in  Figure  3. 

Section  8  describes  the  simplest  possible  mobile  network:  two  links  with  only  one  mobile  node  that 
travels  along  the  x-axis  track.  This  network  determines  baseline  scenarios.  The  simplest  possible  power 
control  simply  assigns  constant  power  to  both  transmitters.  The  best  possible  power  control  assumes  an 
all-knowing  network  “genie”  capability  of  perfect  measurements  and  instantaneous  control.  Good  power 
control  algorithms  will  have  performance  closing  in  on  the  best  possible  power  control  but  must  outper¬ 
form  the  simplest  power  control.  The  mobile  links  were  baselined  at  300  MHz,  1.9  GHz,  and  2.4  GHz. 

The  frequency  band  at  1.9  GHz  was  selected  to  test  the  algorithms.  The  associated  bandwidth  of  1.3  MHz 
implies  the  wireless  link  operate  with  a  frame  duration  of  one  millisecond.  This  frame  rate  is  the  update 
rate  used  for  all  the  wireless  power  algorithms.  That  is,  each  wireless  power  control  algorithm  travels 
along  the  same  mobile  track,  updates  power  every  millisecond,  measures  the  same  noise  power,  is  given 
perfect  interference  measurements,  and  is  subject  to  the  same  system  delays.  These  constraints  enable 
side-by-side  comparisons  of  the  algorithms. 

Section  9  applies  the  fundamental  Foschini  and  Miljanic  (FM)  algorithm  [22]  to  the  mobile  network  at 
1.9  GHz  using  a  1.3-MHz  bandwidth.  Each  link  Z  =  1,  . . . ,  L  updates  its  power  control  command  as 

PC,l(fn+ 1)  =  PC,l(tn) 

The  power  updates  by  the  ratio  of  the  requested  SINR  yx  to  the  measured  SINR  7a/ (in).  If  the  FM  al¬ 
gorithm  converges  (7m (tn)  —>  jx),  the  power  settles  to  a  steady  state  (pcj{tn)  — »  constant).  The 
convergence  rate  is  controlled  by  the  stepsize  f.  The  FM  algorithm  handles  delays  out  to  8  frames  as  the 
mobile  network  traveled  along  its  x-axis  track.  In  comparison  to  the  best  possible  performance  delivered 
by  the  network  “genie,”  the  FM  algorithm  requires  roughly  the  same  power  but  suffers  a  50%  reduction  in 
throughput.  This  loss  of  throughput  is  the  “cost-of-business”  associated  with  a  distributed  algorithm  versus 
the  idealized  centralized  controller.  However,  the  FM  algorithm  breaks  down  when  delay  exceed  8  frames. 
Increasing  delays  cause  power  to  increase  and  throughput  to  decrease. 

Section  10  applies  Multiple  Objective  Distributed  Power  Control  Algorithm  (MODPC)  of  Elumsrati, 
Jantti,  and  Koivo  [20]  to  this  mobile  simulation.  The  MODPC  algorithms  seeks  a  trade-off  between  the 


\  l-d 

IX  \ 

lM,l(tn)  J 


7 


two  objectives  of  power  consumption  and  QoS.  The  MODPC  power  control  is  [20,  Equation  7] 


PC,l{tn+ 1) 


PC,l(tn) 


ApPmin  +  Aq7x 
a PPC,l{tn )  +  A QlM,l{tn )  ’ 


Ap,  Aq  >  0;  Ap  +  Aq  =  1, 


where  the  Ap  «  1  weights  the  trade-off  for  power  savings  while  \p  ~  1  seeks  QoS.  This  solution  contains 
the  FM  algorithm  for  Ap  =  0.  The  MODPC  algorithm,  even  when  weighted  for  QoS,  still  only  handles 
delays  out  to  8  frames.  When  delays  exceed  8  frames,  the  MODPC  algorithm  breaks  down  like  the  FM 
algorithm — power  increases  while  throughput  drops  to  zero. 

Section  1 1  reviews  the  Kalman  filter  as  a  special  case  of  the  H2  algorithms.  Section  12  assesses  two 
II1  algorithms  along  the  same  mobile  network.  Average  power-rate  performance  is  simulated  for  the  mo¬ 
bile  network  at  moderate  delays.  Both  H2  algorithms  breakdown  when  system  delay  exceeds  8  frames. 

The  preceding  sections  show  the  limitations  of  the  FM  algorithm,  Kalman  filtering,  and  the  H°°  al¬ 
gorithms  that  model  time  delay  as  a  noise  source  or  external  disturbance.  Section  13  applies  a  time-delay 
compensation  scheme  to  produce  a  simple  extension  of  the  FM  algorithm.  The  key  idea  is  to  record  the 
power  history  and  estimate  the  “stale”  power  with  an  approximation  from  the  past.  This  simple  extension 
allows  the  FM  algorithm  to  handle  delays  out  to  20  frames. 

Section  14  analyzes  the  effect  of  the  delay  on  the  FM  algorithm.  Even  in  the  simplest  network,  con¬ 
sisting  of  exactly  one  link,  system  delays  cause  instabilities.  Stepsize  selection  enforces  stability  but  only 
handles  a  limited  range  of  delays  and  also  reduces  the  adaptive  speed  of  the  FM  algorithm.  Thus,  the  sig¬ 
nificance  of  the  time-delay  compensation  scheme  is  the  stabilization  of  the  FM  algorithm  without  reducing 
its  adaptive  speed. 

Section  15  assesses  the  simulations.  The  side-by-side  comparison  of  these  distributed  wireless  power 
control  algorithms  operating  in  a  tactical  RF  environment  reveal  similar  performance  across  all  algorithms: 
near-optimal  performance  at  moderate  delays  but  a  breakdown  at  longer  delays,  time-delay  compensation 
stabilizes  the  FM  algorithm.  There  arc  several  points-of-departure  based  these  findings. 

Section  15.1  points  out  that  the  wireless  power  control  algorithms  operated  using  ideal  measurements 
and  synchronized  updates.  Therefore,  assessments  should  encompass  non-ideal  limited  measurements 
available  from  current  radios  and  asynchronous  updates. 

Section  15.2  points  out  that  wireless  power  control  algorithms  only  controlled  power.  Performance  will 
be  enhanced  by  optimizing  power  with  other  control  modes:  beamforming,  channel  selection,  scheduling, 
and  routing.  The  payoff  from  these  extra  control  “knobs”  may  be  substantial  but  exacerbates  system  delay 
and  instability.  The  mobile  wireless  simulations  can  assess  both  the  performance  and  the  stability  under 
the  assumption  of  perfect  cross-layer  communication.  If  the  proposed  algorithms  deliver  significant  per¬ 
formance  shifts  over  the  power-only  control,  consideration  of  the  messaging  complexity  is  the  next  step  to 
assess  the  feasibility  of  these  other  control  modes. 

Section  15.3  points  out  that  the  wireless  power  control  algorithms  operated  over  a  narrowband  channel. 
Section  15.3  presents  a  20- MHz  channel  centered  at  1.9  GHz  from  the  mobile  simulations.  This  channel 
is  frequency  selective.  Although  the  preceding  wireless  power  control  algorithms  can  operate  over  such 
channels,  the  current  approaches  point  to  an  OFDM-CDMA  waveform  to  handle  this  channel.  Regardless 
of  the  waveform,  the  wideband  mobile  wireless  simulations  can  compare  performance  and  the  stability  of 
both  the  narrowband  algorithms  against  the  OFDM-CDMA  approaches  under  ideal  measurements.  If  the 
OFDM-CDMA  approaches  algorithms  deliver  significant  performance  shifts  over  the  narrowband  algo¬ 
rithms,  robustness  of  the  OFDM-CDMA  approaches  is  the  next  step  for  assessment. 


1.7  NOTATION 


The  real  and  complex  numbers  arc  denoted  by  M  and  C,  respectively.  Likewise,  the  corresponding  L- 
dimensional  vector  spaces  are  denoted  RL  and  CL,  respectively.  The  transpose  of  x  £  ML  is 


The  Hermitian  of  x  G  CL  is 


=  [xi, 


^  =  [xi, 


Xl\  • 
Xl\  , 


where  the  overline  denotes  the  complex  conjugate.  A  matrix  A  G  RMxN  then  A  is  said  to  have  size  M  x 
N. 


Table  2.  Network  functions. 


Notation 

Description 

P  r(t) 

P  R.(t) 

P  N(t) 

P  i(t) 

P  cW 
7m  (t) 

7 x(t) 
Rx{t) 

A 

A  tD(t) 

Transmitted  power  vector  (W) 
Received  power  vector  (W) 
Noise  power  vector  (W) 
Interference  power  vector  (W) 
Control  power  vector  (W) 
Measured  SINR  vector 

External  SINR  vector 

External  throughput  vector 
Uplink  delay  (m-sec) 

Downlink  delay  (m-sec) 

Table  3.  RF  functions  and  parameters. 


Notation 

Description 

fc 

Carrier  frequency  (MHz) 

Ib 

Bandwidth  (MHz) 

G 

Channel  gain  matrix 

Gs 

Channel  signal  matrix:  Gs  =  diag(G) 

Gi 

Channel  interference  matrix:  Gj  =  G  -  Gs 

Rb 

Bit  rate  (bps) 

sT(t) 

Transmitted  signal  vector  (V) 

S  n{t) 

Received  signal  vector  (V) 

tF 

Frame  duration  (ms) 

Table  4.  Acronyms. 


Description 

BPSK 

CDMA 

TDMA 

SINR 

Binary  Phase-Shift  Keying 

Code  Division  Multiple  Access 

Time  Division  Multiple  Access 
Signal-to-Noise  and  Interference  Ratio 
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2.  THE  NARROWBAND  NETWORK  MODEL 

This  section  presents  the  signal,  interference,  and  noise  models  for  a  narrowband  wireless  network. 
The  gain  matrix  models  the  signal  power  and  the  network’s  self-interference.  The  noise  power  is  set  by 
the  external  RF  noise.  The  networks’s  signal-to-interference  and  noise  ratio  (SINR)  determines  an  upper 
bound  on  the  throughput.  Thus,  the  network  performance,  as  measured  by  the  power  and  rate,  are  deter¬ 
mined  by  the  gain  matrix  and  the  noise. 

2.1  POWER  FLOWS  IN  THE  NETWORK 

Figure  4  is  a  schematic  of  a  two-link  network.  Each  transmitter  is  assigned  a  designated  receiver  de¬ 
fined  by  the  link  (green  lines).  All  links  share  the  same  frequency  band  and  all  broadcast  at  the  same  time. 
These  simultaneous  broadcasts  constitute  the  interference  (red  lines)  exacerbating  the  additive  in-band 
noise  also  impinging  upon  each  receiver.  For  a  general  network,  let  the  links  be  indexed  l  =  1,  . . . ,  L. 


Figure  4.  Two-link  network  schematic. 


Assume  the  receivers  intercept  the  transmitter’s  signals  according  to  the  following  narrowband  model : 


sr,  i(f) 

_sr,l  (f  )_ 

hn(t) 

/112(f)  • 

. .  hlL{t) 

h2\ (t) 

/122(f)  • 

■■  fl2L(f) 

hLi(t) 

hL2(t )  . 

•  •  hLi{t) 

ST,1  (f ) 

'm  (f)" 

+ 

VL(t)_ 

In  the  notation  of  Table  3,  the  channel  matrix  H{t )  links  the  transmitted  signal  vector  to  the  received  sig¬ 
nal  vector: 

sR(t)  =  H(t)sT(t) +  rj(t).  (1) 

Signal  power  and  interference  power  are  computed  from  Equation  (1).  The  received  power  is 

PR,i(t)  =  E[\sRj(t)\2} 

L  L 

=  e'' 


+  E 


EE  him  (t)ST.m  (t)him'  (f  ) 

jn=  1  m'= 1 
L 

E  (him(t)sT,m{t)vi(t))  +E[\rn(t)\ 2 


,171=1 
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Assume  the  channel  is  independent  of  the  signals  and  the  noise.  Assume  the  noise  is  independent  of  the 
signals  and  is  zero  mean.  Then  the  cross  terms  vanish  and  the  received  power  simplifies  as 

L  L 

PR,l(t )  =  EEe  [^Zra  {t)him>(t)]  E 

m=l  m'= 1 

If  the  signals  arc  uncorrelated,  their  cross  terms  vanish  and  the  received  power  exposes  the  signal,  interfer¬ 
ence,  and  noise  as  follows: 


■st  ,m  ( t)sT  ,mf  (*)  +PN,l(t). 


PR, =  ^2  E  [\hlm(t)\2]  E  [\sT,m{t)\2]  +  PN,l{t) 


m=  1 
L 


^2  9lm{t)pr,m{t )  +PN,l(t) 


m=  1 


signal 


interference 


=  9ll(t)pr,l(t )  +  £<U*)Pr,m(f)  +  PN,l(t )• 

m^l 


PR, lit) 

9u{t)  ■ 

••  9lL{t) 

PT,l{t ) 

PN,l(t) 

p Bit)  = 

= 

+ 

JPR, Lit) _ 

9Ll(t)  . 

■■  9LI,{t)_ 

Pt,l  {t)_ 

PN,L(t)_ 

(2) 


Consequently,  the  network  performance  is  governed  by  gain  matrix  G  and  the  noise  level: 

=  G(t)pT(t)  +  PjvW-  (2) 


Equation  (3)  is  the  basis  for  the  wireless  power  control  algorithms.  This  equation  is  based  on  the  following 
assumptions: 


Assumptions  1  Narrowband  Network  Model 

NB-1  Links  transmit  simultaneously  in  the  same  frequency  band  (no  TDMA) 
NB-2  Link  assignments  do  not  change  (no  handoffs) 

NB-3  Channels  are  narrowband  (multiplicative  gain  matrix) 

NB-4  Channels  are  independent  of  the  signals 

NB-5  Channels  are  independent  of  the  additive  link  noise 

NB-6  Noise  is  zero  mean,  Gaussian,  and  spectrally  flat 

NB-7  Noise  is  uncorrelated  with  the  signals 

NB-8  Signals  are  uncorrelated  between  the  links 

NB-9  Signals  “occupy”  the  entire  band  with  approximately  flat  spectra. 


These  assumptions  do  not  include  the  time  delays  associated  with  the  measurements  and  signal  pro¬ 
cessing  necessary  to  support  the  wireless  power  control  systems.  The  generalization  of  Equation  (3)  to  in¬ 
clude  time  delays  is  discussed  in  Section  5.  The  signal  power,  the  interference  power,  and  the  noise  power 
determine  the  SINR  The  SINR  sets  an  upper  bound  on  link  capacity  and  is  reviewed  in  Section  2.3.  The 
noise  is  modeled  in  the  next  subsection. 
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2.2  NOISE  POWERS  AT  THE  RECEIVER 


The  noise  level  in  the  receiver  should  be  set  at  the  detector  of  the  receiver.  A  standard  noise  model  at 
the  input  of  the  detector  assumes  the  noise  is  the  sum  of  the  external  RF  noise  observed  through  the  re¬ 
ceiver  and  the  receiver’s  internal  noise.  Figure  5  illustrates  the  power  flows  arriving  at  the  receiver’s  an¬ 
tenna  and  into  its  detector.  The  RF  signal,  the  external  RF  noise,  and  the  interference  all  travel  through  the 
receiver’s  RF  chain.  This  report  assumes  the  external  RF  noise  dominates  the  receiver’s  noise.  Therefore, 
the  external  SINR  at  the  antenna  input  is  the  same  as  the  SINR  at  the  input  to  the  detector. 


Figure  5.  Power  flows  into  the  receiver:  pt,rf  transmitted  signal  power  arriving  at  the  antenna;  pn,rf 
environmental  RF  noise  power;  pi,rf  interference  RF  power;  pn,lo  internal  or  local  noise  power. 


The  external  RF  noise  is  modeled  as  zero-mean,  complex-valued  Gaussian  with  flat  spectrum  across 
the  frequency  band  (see  NB-6).  As  such,  only  the  variance  or,  equivalently,  in-band  power  level  pn,rf,  is 
required  to  model  this  Gaussian  noise.  This  noise  level  is  modeled  as  [39,  page  76] 

PN,  RF  =  KsTeqfB,  (4) 

where  Boltzmann’s  constant  is  kb  =  1.38  x  10-23  (W  •  K_1  •  Hz' 1 ),  Teq  is  the  effective  operating 
temperature  (K),  and  fs  is  the  bandwidth  (Hz).  Figure  6  plots  Teq  as  a  function  of  frequency  for  selected 
RF  noise  sources.  Equation  (4)  models  the  narrowband  noise  or,  equivalently,  assumes  a  flat  spectrum. 


Figure  6.  Sky  temperatures  Teq  (Kelvins)  generated  by  background  noise  sources  [39,  page  74], 
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Wideband  noise  is  modeled  by  integrating  the  noise  spectral  density  [39,  page  76]: 


PnU)  =  KfiTeq(/)  [watts  •  Hz  x] 


over  the  frequency  band 


Pn 


rfc+fB/z 
’  Sc-fs!^ 


PN(f)df 


[dBW]. 


Figure  7  plots  the  urban  noise  spectral  density  centered  at  fc  =  300  MHz  over  the  frequency  band 
fs  =  10  MHz.  The  plot  shows  near-constant  noise  density  generates  in-band  noise  power  of  ~  —106 
[dBW].  Therefore,  the  flat-spectrum  of  Assumption  NB-6  is  credible  for  relative  bandwidths 


Sb 

fc 


10 

300 


«3%. 


For  convenience,  the  noise  levels  used  in  the  forthcoming  network  simulations  are  listed  in  Table  5. 


Noise  Power:  Urban  Man-Made 


Figure  7.  Noise  spectral  density  of  the  10-MHz  band  centered  at  300  MHz. 


Table  5.  Urban  and  suburban  noise  levels  (dBW)  for  selected  bandwidths  fB  and  carrier  frequencies  fe¬ 


Frequency 

Bandwidth 

Urban 

Power 

Suburban 

Power 

300  MHz 

25  kHz 

-132 

-147 

300  MHz 

1  MHz 

-116 

-131 

300  MHz 

10  MHz 

-106 

-121 

1.9  GHz 

1.3  MHz 

-133 

-149 

1.9  GHz 

10  MHz 

-124 

-140 

2.4  GHz 

10  MHz 

-126 

-142 

2.4  GHz 

20  MHz 

-123 

-139 

ll 


2.3  SINR  AND  THROUGHPUT 


Assumptions  1  model  the  received  power  on  the  Z-th  link  by  Equation  (2)  and  repeated  here  for  conve- 


interference 


PR, lit)  =  gu(t)pr,l(t)  +  ^9lm(t)PT,m(t)  +  PN,l(t)  . 

m^l 


Definition  1  (SINR)  Under  Assumptions  1,  the  signal-to-interference-and-noise  ratio  (SINR)  measured  on 
the  l-th  link  is 

._  9u(t)pr,i(t ) 


7 M,i{t)  ■= 


PN,l(t )  +  Y,m#9lm(t)pT,m(ty 


The  signal-to-interference-and-noise  ratio  (SINR)  governs  the  throughput  for  specific  waveforms  and  also 
the  maximum  capacity  or  Shannon  capacity  that  a  link  can  support.  This  report  uses  the  Shannon  capacity 
to  bound  best  possible  performance  shifts  that  wireless  power  control  algorithms  deliver  rather  than  getting 
entangled  with  specific  waveforms,  detection  schemes,  estimation  issues,  and  hardware  limitations.  The 
Shannon  capacity  for  a  single  link  that  is  ruled  by  the  signal-to-noise  ratio  (SNR)  [27,  Equation  4.1] 

Rb  =  fB  log2(l  +  SNR), 

provided  the  noise  is  zero-mean,  flat  spectrum,  Gaussian  that  is  uncorrelated  with  the  signal.  If  the  inter¬ 
fering  signals  “look  like”  this  flat-spectrum  Gaussian  noise  and  arc  also  uncorrelated  with  the  link’s  signal, 
the  capacity  generalizes  from  the  SNR  to  the  SINR. 


Definition  2  (Narrowband  Throughput)  Under  Assumptions  1,  the  throughput  on  the  l-th  link  is  the 
Shannon  capacity  integrated  across  the  common  bandwidth 

Rb,l{t )  :=  fB  log2(l  +  7 M,lit)). 


Thus,  the  network  performance — power  and  throughput — is  ruled  by  the  gain  matrix  G(t)  and  the  noise 
power  p pf(t).  The  diagonal  of  the  gain  matrix  determines  the  signal  power.  The  off-diagonal  terms  of  the 
gain  matrix  determine  the  interference.  Networks  using  Code-Division  Multiple  Access  (CDMA)  suppress 
interference.  Roughly  speaking,  CDMA  “coding  gain”  divides  the  interference.  The  next  section  details 
this  coding  gain  for  use  in  the  wireless  simulations. 


2.4  CDMA  CODING  GAIN  AND  THROUGHPUT 

The  gain  matrix  controls  the  signal  and  interference.  Coding  gain  suppresses  the  interference  and 
boosts  the  SINR  and  throughput.  Code-Division  Multiple  Access  (CDMA)  is  one  approach  using  signal¬ 
processing  to  boost  the  SINR  [52]:  “In  a  CDMA  system,  all  users  transmit  information  over  the  same  time 
and  frequency,  and  each  user  signal  is  recovered  through  a  particular  code  sequence.”  The  received  signal 
on  l-th  link  adapts  from  [6,  Equation  2]  as  follows: 

L 

sr, i  (Pt!  t)  =  £  hitky/PT,khsi-  (5) 

k= 1 
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Under  Assumptions  1,  the  channel  matrix  H  determines  the  gain  matrix  G  =  E\\H\2}  and  includes  the 
physical  RF  channel,  the  transmit  antenna,  the  receive  antenna,  and  the  processing  losses.  For  Binary 
Phase-Shift  Keying  (BPSK),  the  information  bit  is  bk  €  {  —  1,1}.  The  symbol  vector  s/,.  typically  derives 
from  a  pseudo-noise  sequence.  For  example,  if  the  receiver  on  the  Z-th  link  uses  a  decision  vector  d/,  the 
SINR  adapts  from  [6,  Equation  7]  as  follows: 

9ll(^hsl)PT,l 


7KPt) = 


PN,l\\Al\\2  +  9lk(dl,Sk)pT,k 


If  d/  =  S/,thc  SINR  simplifies  to 


7/(Pt)  = 


9llPT,l 


PN,l  +  9lk(Sl,Sk)pr,k 

Denote  the  lower  bound 

9c  :=  min{|(s/,sfc)|  :  k,l  G  {1....L}}. 

The  SINR  admits  the  bound  [27,  Equation  14.5]: 

/  ,  ,  9llPT,l 

7z(Pt)  < 

In  terms  of  the  gain  matrix 

G  =  GS  +  Gr,  Gs  = 


PN,l  +  5c1  9lkPT,k 


9n  o 
0  922 


0 

0 


Gr  =  G  -  Gs, 


0  0  ...  gLL . 

coding  gain  reduces  the  interference  by  changing  the  gain  matrix: 

Gt-GS  +  -G;. 

9C 

Simulations  that  list  “coding  gain”  imply  that  Equation  (6)  bounds  the  SINR. 


(6) 


For  the  throughput,  the  capacity  equation  must  be  specialized  to  BPSK-CDMA  coding.  Assume  each 
network  transmits  CDMA  symbol  vectors  consisting  of  Nc  chips  of  duration  Tc.  The  network  bandwidth  is 


I'b  = 


1 

Tc’ 


each  symbol  vector  has  duration  Tsym  =  Nc  x  Tc  ,  and  the  symbol  rate  is 


j1  N  x  T 

^-sym  iVc  *  -Lc 


Nc~ 


The  throughput  specialized  for  BPSK-CDMA  is  adapted  from  [6,  Equation  5].  Each  symbol  vector  carries 
one  information  bit  by  Equation  (5): 


Rb  —  Rsym  x  (1  ex]> (  7m)) 


x  (1  -  exp(-7 M))Nc  [bps]. 


(V) 


This  adaptation  assumes  no  overhead  symbols  (L  =  1)  and  the  standard  efficiency  function  [6].  This 
BPSK-CDMA  throughput  is  dominated  by  the  standard  capacity  equation.  Therefore,  whenever  coding 
gain  is  invoked,  throughput  can  be  computed  using  the  BPSK-CDMA  rate  of  Equation  (7)  or  throughput 
can  be  bounded  by  the  standard  capacity  equation.  The  simulations  that  invoke  the  coding  gain  bound  of 
shown  in  Equation  (6)  also  bound  the  throughput  by  standard  capacity  equation — not  the  BPSK-CDMA 
rate  of  Equation  (7). 
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3.  RF  PROPAGATION  IN  URBAN  SCENARIOS 

Figure  8  shows  the  three-dimensional  urban  environment  where  the  mobile  wireless  networks  oper¬ 
ate.  Figure  9  is  the  top  view  of  a  two-link  network  operating  in  this  urban  area.  Both  transmitters  and  one 
receiver  are  at  street  level.  The  second  receiver  is  located  on  the  rooftop  80  meters  above  the  street.  The 
transmitter  on  Link  2  will  travel  along  the  x-axis.  The  view  shows  transmitter  TX-2  at  the  origin.  This 
two-link  network  is  the  smallest  and  simplest  mobile  wireless  network  simulation.  This  simulation  is  the 
common  operating  environment  to  assess  the  forthcoming  wireless  power  control  algorithms. 


Rosslyn,  Virginia 


r  (meters) 


Figure  8.  Urban  propagation  environment. 
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Figure  9.  Top  view  of  the  two-link  network. 

RF  propagation  in  this  simulated  3-D  environment  is  governed  by  the  wave  equation  [24,  Chapter  3]. 
Approximate  solutions  to  the  wave  equation  are  computed  by  Numerical  Electromagnetics  Code-Basic 
Scattering  Code  (NEC-BSC)  [55].  NEC-BSC  approximates  the  electric  and  magnetic  field  solutions  to  the 
3-D  wave  equation  by  reflection  and  diffraction  models.  The  approximation  for  the  electric  field  generated 
by  a  source  sums  the  rays’  interaction  with  the  3-D  environment  to  arrive  at  a  position  r: 

E(r)  =  Einc  +  EP(<E)  +  E  E6(^)  +  E  EW(<JU) 

r&V{  r)  eG£(r)  i>£V(r) 

+higher-order  terms. 

The  first  term  is  the  incident  electric  field  Ejnc  produced  by  the  direct  path  from  the  source  to  posi¬ 
tion  r,  provided  there  is  no  blockage.  The  second  term  sums  each  electric  field  generated  by  the  ray  trav¬ 
eling  from  the  source,  reflecting  off  a  plate  at  position  qr,  and  arriving  at  position  r.  Each  plate  reflection 
models  the  dielectric  of  the  plate  and  consequent  spreading  factor  after  reflection.  The  third  term  sums 
each  electric  field  generated  by  the  ray  traveling  from  the  source,  diffracting  on  an  edge  at  position  qe,  and 
arriving  at  position  r.  Figure  10  shows  the  diffraction  cone  geometry  for  an  edge.  Each  edge  diffraction 
models  the  dielectric  of  the  edge  and  consequent  spreading  factor  after  diffraction.  The  edge  diffraction 
terms  correct  for  the  discontinuities  in  the  incident  and  plate  reflected  fields.  The  fourth  term  sums  each 
electric  field  caused  by  rays  traveling  from  the  source,  diffracting  at  vertex  q,.,  and  arriving  position  r. 
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Each  edge  diffraction  models  the  dielectric  of  the  vertex  and  consequent  spreading  factor  after  diffraction. 
The  vertex  diffraction  terms  correct  for  the  discontinuities  in  the  edge-diffracted  fields. 


Figure  10.  Edge  diffraction  [56,  Figure  2], 

NEC-BSC  differs  from  other  ray  tracing  programs  by  modeling  not  only  the  plate  reflections  but  also 
including  the  diffraction  terms.  The  diffraction  modeling  in  NEC-BSC,  even  at  the  ray  level,  also  differs 
from  other  ray  tracing  programs.  Diffraction  coefficients  arc  derived  from  the  high-frequency  asymptotic 
solutions  of  canonical  problems.  These  asymptotic  approximations  are  more  accurate  as  the  associated 
“large  parameters”  increases.  The  “large  parameters”  arc  generally  a  function  of  the  distances  from  the 
source  to  the  diffraction  point  and  from  the  diffraction  point  to  the  position  r.  Multiple  interactions  be¬ 
tween  flat  faces,  edges  and  vertices  are  needed  for  complex  environments.  Cascading  these  high-frequency 
asymptotic  solutions  determine  these  higher-order  interactions.  Complex  environments  have  an  infinite 
number  of  higher-order  interactions.  Higher-order  terms  compensate  discontinuities  in  first-order  and  other 
lower-order  terms.  Table  6  lists  the  second-order  terms. 

Table  6.  Second-order  diffraction  terms. 


Plate-to-plate 

Plate-to-edge 

Plate-to-vertex 


Edge-to-plate 

Edge-to-edge 

Edge-to-vertex 


Vertex-to-plate 

Vertex-to-edge 

Vertex-to-vertex 


Third-order  terms  follow  a  similar  pattern  (e.g.,  Plate-to-plate-to-plate,  Plate-to-plate-to-edge,  . . . ). 

A  correct  application  of  NEC-BSC  requires  that  the  dominant  terms  in  an  environment  be  determined  to 
guarantee  accuracy  of  the  solution.  Such  a  convergence  study  was  undertaken  to  support  the  propagation  in 
[63].  Table  7  lists  the  diffraction  terms  necessary  to  approximate  the  electric  field  to  30-dB  accuracy.  The 
dominant  terms  were  limited  to  plate  and  edge  interactions. 
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Table  7.  Diffraction  terms  for  the  urban  simulation. 


Order 

Terms 

0 

Incident 

1 

Plate,  edge 

2 

Plate-to-plate,  Edge-to-plate,  Edge-to-edge 

The  reflections  and  diffractions  require  the  dielectrics  of  the  building  and  the  ground.  Table  8  describes 
the  modeling.  The  buildings  are  modeled  with  roofs  and  sides  as  concrete  slabs  0.3  meters  thick.  Win¬ 
dows,  iron  reinforcements,  or  rooftop  materials  are  not  modeled.  The  ground  plane  is  modeled  as  a  hal- 
plane  of  wet  earth.  Variations  in  elevation  arc  not  modeled.  The  wet-earth  dielectric  was  selected  to  match 
the  modeling  in  Daly  [15].  Wet  earth  absorbs  more  RF  energy  than  the  asphalt  listed  for  comparison. 

Table  8.  Urban  dielectrics. 


Slab 

thickness 

Relative 

Permittivity 

Conductivity 

References 

Buildings 

0.3  (m) 

6.0 

0.01 

[15] 

Wet  earth 

Half  space 

25 

0.005 

[15] 

Asphalt 

Half  space 

2.56 

0.0011 

[65,  4-30] 

This  simple  3-D  urban  modeling  produces  representative  RF  environments.  These  RF  environments 
model  large-scale  multipath  behaviors  rather  than  absolute  levels  [63].  As  such,  these  simulations  arc  best 
employed  to  reveal  relative  performance  shifts  between  the  wireless  power  control  algorithms. 

Performance  of  any  wireless  power  control  algorithm  is  ruled  by  the  gain  matrix.  Section  3.1  starts 
the  discussion  by  introducing  the  channels  generated  by  the  Transmitter  TX-1  broadcasting  at  300  MHz 
(Figure  9).  Section  3.2  shows  how  the  channel  plots  vary  as  Transmitter  TX-2  travels  through  the  urban 
environment  along  the  x-axis  while  broadcasting  at  1 .9  GHz.  The  channel  plots  arc  then  combined  to  pro¬ 
duce  the  gain  matrix  for  the  wireless  network  generated  by  the  mobile  wireless  network  with  TX-2  travel¬ 
ing  along  the  x-axis  track.  Section  3.3  provides  the  associated  gain  matrix  at  2.4  GHz.  Although  the  3-D 
urban  model  is  simple,  the  gain  matrices  show  large-scale  variations  caused  by  the  street  geometry  (large- 
scale  shadowing)  and  fine-scale  variation  (fast  fading)  caused  by  the  local  diffraction  effects.  Any  wireless 
power  control  algorithm  must  adapt  to  both  types  of  fading. 
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3.1  URBAN  PROPAGATION  AT  300  MHZ 


Figure  1 1  is  a  top  view  looking  down  on  the  street-level  propagation  produced  by  transmitter  TX- 1 
broadcasting  at  300  MHz  or  with  a  wavelength  of  1  meter.  The  color  map  is  the  channel  gain  from  TX-1  to 
street-level  receivers  produced  by  a  NEC-BSC  simulation.  Both  the  transmitter  and  receiver  operate  with 
a  half-wavelength  vertical  dipole.  Both  the  transmitter  and  receiver  are  “street  level”  or  that  their  antennas 
are  placed  2  meters  above  the  ground  plane  (z  =  2). 


Channel  at  /=  0.3  GHz;  RX  z-2  (m) 
Average— 109  dB;  Max=-18  dB 


x  (m):  Buildings:  s=  6,  tan(5)=0. 099909;  Wet  Earth:  s=  25,  tan(5)=0. 0047956 


Figure  1 1 .  Top  view  of  the  300-MHz  channel  from  transmitter  TX-1 . 

Even  in  this  simple  3-D  model,  the  channel  gain  shows  detailed  interference  patterns  near  the  trans¬ 
mitter.  The  plot  also  shows  RF  energy  diffracting  around  the  buildings  and  the  resulting  shadow  zones 
caused  by  each  building.  In  these  low-power  shadow  zones,  the  discontinuities  caused  by  using  only  the 
second-order  model  (Table  8)  arc  visible.  The  low-power  fields  could  be  smoothed  by  adding  higher-order 
terms — incurring  a  massive  computational  expense — to  approximate  the  gain  down  to  50  dB  rather  than  30 
dB. 
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The  next  three  figures  display  the  effect  of  placing  a  transmitter  or,  equivalently,  a  receiver  to  a  rooftop. 
Figure  12  shows  the  3-D  urban  environment  with  a  transmitter  on  the  rooftop  that  is  80  meters  above  the 
street.  The  transmitter’s  antenna  is  2  meters  above  the  rooftop.  Figure  13  shows  the  rays  computed  by 
NEC-BSC  that  link  the  rooftop  transmitter  to  a  street-level  receiver.  Close  examination  shows  that  most 
of  the  ray  paths  include  a  diffraction  off  a  building  edge. 

Figure  14  is  a  top  view  looking  down  on  the  street-level  propagation  produced  by  the  rooftop  transmit¬ 
ter  broadcasting  at  300  MHz  down  to  street-level  receivers.  The  transmitter  is  located  at  x  =  75,  y  =  50, 
and  z  =  82  meters.  Both  transmitters  and  receivers  employ  a  half-wavelength  vertical  dipole.  The  color 
map  uses  a  different  scale  than  Figure  1 1  to  better  show  the  relative  uniform  coverage,  or  equivalently,  the 
tactical  communication  advantage  that  this  location  offers. 

Although  the  channels  in  Figure  14  are  computed  using  the  rooftop  transmitter  to  a  street-level  re¬ 
ceiver,  reciprocity  implies  this  plot  also  shows  the  channels  from  street-level  transmitters  to  the  rooftop 
receiver  RX-2  of  the  two-link  network  shown  in  Figure  9.  In  particular,  this  simulation  models  street-level 
transmitters  traveling  through  this  environment  whose  signal  is  intercepted  by  RX-2.  The  next  section  ex¬ 
amines  specific  channels  along  the  track  followed  by  TX-2. 


Rosslyn,  Virginia 


400 


Figure  12.  Rooftop  transmitter  in  an  urban  propagation  environment  (transmitter’s  height  exaggerated  for 
clarity). 
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10  dB 


Figure  13.  Rays  emanating  from  the  rooftop  transmitter  down  to  a  street-level  receiver. 


Channel  at  /=  0.3  GHz;  RX  2=2  (m) 
Average=-1 1 7  dB;  Max=-63  dB 
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Figure  14.  Top  view  of  the  300-MFIz  channel  from  the  rooftop  transmitter  (x  =  75,  y  =  50,  z  =  82  meters) 
down  to  street-level  receivers. 
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3.2  URBAN  PROPAGATION  AT  1.9  GHZ 


Figure  15  shows  the  smallest  and  simplest  possible  mobile  wireless  network  consisting  of  four  nodes 
or  two  links  and  only  one  mobile  node.  The  mobile  node  carries  transmitter  TX-2  along  the  x-axis  track. 
All  nodes  employ  half-wavelength  vertical  dipole  antennas. 


Rosslyn,  Virginia 


Figure  1 5.  Top  view  of  the  mobile  network;  TX-2  travels  along  the  x-axis. 

The  next  four  plots  show  the  channel  generated  by  TX-2  at  selected  locations  along  its  x-axis  track: 
x  =  0,  50,  100,  and  150  meters.  The  plots  show  the  power  radiating  from  TX-2  flowing  through  the  city  up 
to  second-term  accuracy.  Consequently,  the  high-gain  regions  are  a  better  approximation  to  the  gain  than 
the  low-gain  regions  in  the  narrow  side  streets. 
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Channel  at  /=  1 .9  GHz;  RX  z=2  (m) 
Average=-129  dB;  Max=-34  dB 


x  (m):  Buildings:  s=  6,  tan(8)=0.015775;  Wet  Earth:  G=  25,  tan(S)=0.0007572 

Figure  1 6.  Coverage  of  TX-2  at  x  =  0  meters. 


Channel  at  /=  1 .9  GHz;  RX  z=2  (m) 
Average=-1 19  dB;  Max=-34  dB 


x  (m):  Buildings:  s=  6,  tan(8)=0.015775;  Wet  Earth:  s=  25,  tan(8)=0.0007572 


Figure  17.  Coverage  of  TX-2  at  x  =  50  meters. 
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\h\l  dB  \h\2  dB 


Channel  at  /=  1 .9  GHz;  RX  z=2  (m) 
Average=-1 19  dB;  Max=-33  dB 


x  (m):  Buildings:  S=  6,  tan(8)=0.015775;  Wet  Earth:  s=  25,  tan(S)=0.0007572 

Figure  18.  Coverage  of  TX-2  at  x  =  100  meters. 


Channel  at  /=  1 .9  GHz;  RX  z=2  (m) 
Average=-1 17  dB;  Max=-34  dB 


x  (m):  Buildings:  e=  6,  tan(8)=0.015775;  Wet  Earth:  s=  25,  tan(S)=0.0007572 


Figure  19.  Coverage  of  TX-2  at  x  =  150  meters. 
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\h\l  dB  \h\2  dB 


If  the  channels  of  this  two-link  mobile  network  are  collected  as  TX-2  travels  along  the  x  axis,  the  re¬ 
sulting  gain  matrix  has  the  form 


TX— 1  TX-2 


G{x)  = 


RX-1 

911  012(2;) 

'-102 

012  (x) 

RX-2 

021  022 (x) 

-93 

022  (x)_ 

[dB], 


The  first  column  is  constant  because  transmitter  TX- 1  and  both  receivers  are  stationary.  Figure  20  plots  the 
second  column  of  this  gain  matrix  showing  the  link  <722  (x)  and  interference  <712  (x)  terms. 


Rosslyn  at  1900  MHz:  Moving  TX-2;  Gain  Matrix 
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Figure  20.  Gain  matrix  at  1 .9  GHz  for  the  two-link  mobile  wireless  network. 


The  channel  gain  522  (x)  models  the  link  from  the  moving  street-level  transmitter  TX-2  to  its  rooftop 
receiver  RX-2.  This  channel  is  relatively  flat  because  there  is  little  blocking  by  buildings  along  the  track 
(Figure  14).  However,  the  channel  gain  <712  (x)  models  the  interference  from  at  receiver  RX-1  produced  by 
the  moving  street-level  transmitter  TX-2.  The  preceding  plots  show  that  TX-2  broadcasts  up  the  street  con¬ 
taining  RX-1  in  the  initial  part  of  its  track  (Figures  16  and  17).  Once  the  transmitter  is  past  the  entrance  to 
this  street,  the  buildings  block  the  interference  (Figures  18  and  19).  Figure  20  shows  the  large-scale  shad¬ 
owing  caused  by  the  street  geometry  and  the  fast  fading  caused  by  the  local  diffraction  effects. 
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Figure  20  shows  the  channel  h 22(2;)  has  two  distinct  fading  regimes: 


•  x  <  100  meters:  slow  fading;  characteristic  of  two-path  propagation 

•  150  <  x  meters:  fast  fading;  no  strong  path. 

Slow  and  fast  fading  affect  wireless  power  control  algorithms  that  allocate  power  based  on  “predicted”  fad¬ 
ing.  Consequently,  channels  with  longer  spatial  correlation  arc  more  amenable  to  prediction.  The  channel’s 
spatial  correlation  function  is  defined  as  [4,  page  12] 

Chh( Ax)  =  i  Rhh{ Ax)  =  E[h22{x)h22(x  -  Ax)]. 

The  channel’s  correlation  distance  A x(p)  is  defined  as  the  first  positive  distance  where  correlation  drops 
below  the  specified  correlation: 

A x(p)  =  argrnin  {Ax  >  0  :  Chh{ Ax)  <  p}  . 

Figures  21  and  22  compare  the  estimated  spatial  correlation  function  for  slow  and  fast  fading,  respectively 
(7/22 (X)  in  the  fading  estimate  means  that  the  linear  trend  was  removed  from  the  channel  magnitude  be¬ 
fore  estimating  the  correlation).  The  spatial  correlation  function  over  the  slow-fading  region  shows  typical 
“two-path”  fading  pointing  to  Rician  fading.  The  slow-fading  channel  retains  90%  correlation  for  almost 
a  meter:  Ax(90%)  =  0.9  meters.  The  fast-fading  region  rapidly  decorrelates  and  has  little  recognizable 
sidelobe  structure:  Ax(90%)  =  0.07  meters.  Table  9  compares  the  spatial  correlation  distances  for  both 
regions. 


Table  9.  Spatial  correlation  distance:  1.9  GHz. 


Region 

Ax(50%)  Ax(90%) 

Slow 

Fast 

2.32  m  0.90  m 

0.36  m  0.07  m 

Figures  23  and  24  estimate  the  fading  distribution  for  slow  and  fast  fading,  respectively.  The  slow  fad¬ 
ing  registers  Rician  fading  but  identifying  fast  fading  is  more  problematic.  The  lognormal  distribution  is 
best  matched  for  wireless  power  control  based  on  the  Kalman  filter.  Wireless  power  control  algorithms  that 
adapt  to  multiple  fading  regimes  promise  greater  performance  than  algorithms  restricted  to  a  single-mode 
fading. 
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correlation  correlation 


Figure  21.  Spatial  correlation  of  h22(x)  on  -20  <  x  <  100  (m). 


Figure  22.  Spatial  correlation  of  h22(x)  on  150  <  x  (m). 
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Rosslyn  1.9  GHz:  GoF g,,(*)  x  <  100  meters 
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Figure  23.  Goodness-of-Fit  for  \h22(x)\2  on  -20  <  x  <  100  (m). 
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Figure  24.  Goodness-of-Fit  for  \h22(x)\i 2  on  150  <  x  (m). 
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3.3  URBAN  PROPAGATION  AT  2.4  GHZ 


Figure  25  shows  the  gain  matrix  at  2.4  GHz  of  the  two-link  mobile  network  as  TX-2  travels  along  the 
.x-axis.  This  plot  is  included  for  comparison  to  the  1.9-GHz  gain  matrix.  As  in  the  300-MHz  simulation 
and  the  1.9-GHz  simulation,  both  transmitters  and  receivers  employ  a  half-wavelength  vertical  dipoles. 
Both  plots  show  the  same  structure:  large-scale  shadowing  caused  by  the  street  geometry  and  fast  fading 
caused  by  the  local  diffraction  effects.  Any  adaptive  wireless  power  control  algorithm  must  handle  both 
types  of  fading. 


Rosslyn  at  2400  MHz:  Moving  TX-2;  Gain  Matrix 


0  50  100  150  200 

x  (meters):  g  =-105  dB:  g  =-91  dB 

Figure  25.  Gain  matrix  at  2.4  GHz  for  the  two-link  mobile  wireless  network. 
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4.  INTERLUDE:  TRADE-OFFS  AND  PARETO  POINTS 


This  section  reviews  a  few  results  multi-objective  optimization  necessary  to  state  optimal  trade-offs  be¬ 
tween  competing  objective.  Trade-offs  between  competing  objectives  are  conveniently  stated  using  vector 
inequalities. 

Definition  3  The  vector  x  €  R;V,  is  said  to  be  element-wise  non-negative  if  no  element  ofx  is  negative: 
Xm  >  0  denoted  as  x>0.  If  y  £  MA',  the  vector  inequality  x<y  means  y  —  x>0. 

The  set  of  all  element-wise  nonegative  vectors  is  denoted  as 

,=  jy  £  ;  y>0}  . 

This  closed,  convex  orthant  of  R;V  determines  the  geometry  of  this  vector  inequality. 

Definition  4  (Minimum)  [3,  page  45]  Let  Y  C  R;V.  A  point  ymin  £  V  is  called  the  minimum  of  Y  if  and 
only  ifymm<yfor  all  y  £  Y. 

Equivalently,  Y  is  completely  contained  in  the  positive  orthant  R;V  located  at  the  minimum  [3,  page  45]: 

y  c  ymin  +  • 

Figure  26  shows  the  relation  between  a  cone  Y  C  M2  to  the  positive  orthant  R+  for  the  vertex  of  the  cone 
to  be  the  minimum  of  the  cone.  However,  sets  with  a  minimum  are  relatively  rare.  More  common  are  the 
minimum  elements. 


K 


*  min 


Figure  26.  Vertex  of  a  cone. 


Definition  5  (Minimal  Elements)  A  point  yME  £  Y  is  called  a  minimal  element  if  and  only  if  for  all 
y  e  Y  such  that  y<yME  =>  y  =  yME- 

The  geometry  for  yME  £  Y  to  be  a  minimal  element  requires  that  the  negative  orthant  located  at  yME 
intersects  Y  only  at  the  minimal  element.  This  geometry  is  formalized  in  the  Contact  Theorem  [77,  Con¬ 
tact  Theorem].  Figure  27  shows  how  the  negative  orthant  — must  attach  to  Y  c  R2  for  yME  to  be  a 
minimal  element. 
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Theorem  1  (Contact  Theorem  [77])  :  A  point  yME  G  Y  tis  a  a  minimal  element  if  and  only  if 


{yMEl  —  (yME  —  K+) n  Y- 


Figure  27.  A  minimal  element. 

Optimal  trade-offs  for  a  multiobjective  function  h  :  X  C  Mm  ->  rn  arc  determined  by  minimal 
elements  of  the  performance  image.  Figure  28  illustrates  the  image  of  a  Pareto  optimal  point  for  h  :  X  C 

Rm  M2. 


Definition  6  (Pareto  Optimal)  [ 3,  page  ]  Let  h  :  X  C  RM  —?  RN  and 

h(X)  :=  (h(x)  :  x  €  RM}  . 

A  element  xp  G  X  is  called  Pareto  optimal  if  and  only  if  h(xp)  is  a  minimal  element  of  h(X). 


< 


R" 

+ 


y 


Figure  28.  Pareto  optimal  point  for  the  mapping  h:ic  RM  ->•  ]R2. 
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The  image  of  all  Pareto  optimal  points  is  called  the  Pareto  front  or  Pareto  frontier.  Figure  29  illus¬ 
trates  such  Pareto  front  marked  by  the  heavy  line.  The  dashed  lines  indicate  that  this  Pareto  front  does  not 
“wrap”  past  the  vertical  and  horizontal  axis. 


_ {  ^ 

I  ^ 

I 

i 

Figure  29.  Pareto  front  of  the  mapping  h  :  A'  c  RM  — >  R2. 


This  section  closes  with  a  few  observations  about  the  Pareto  front  and  its  specialization  to  the  power- 
rate  image.  The  first  observation  points  out  that  the  Pareto  front  must  lie  in  the  boundary  of  the  perfor¬ 
mance  image 

Lemma  1  Let  h  :  X  C  Mm  — >  R;V.  The  Pareto  front  ofh(X )  lies  in  the  boundary  ofh(X). 


If  the  image  is  its  own  boundary,  such  as  the  graph  of  a  function,  the  analytic  properties  of  the  function 
may  force  the  image  to  be  its  own  Pareto  front. 


Lemma  2  Let  h  :  X  C  M  — y  M  be  continuous  and  strictly  decreasing.  Define  H  :  X  — >  M2  as 


H(x) 


x 

h{x) 


Then  the  Pareto  front  ofH(X)  is  H  (X ). 


When  the  generic  Pareto  theory  is  specialized  to  the  power-rate  image,  the  negative  orthant  is  replaced 
by  the  quadrant  pointing  to  “less  power  and  more  rate”.  Let  QcK2  denote  this  quadrant: 


Q  := 


rz 


■  PT.XT.  >  0 


Figure  30  shows  this  quadrant  in  contact  with  the  Pareto  front  of  the  power-rate  image.  The  upper  bound¬ 
ary  of  the  power-rate  image  is  the  Pareto  front — no  variation  in  the  power  vector  exists  that  simultaneously 
decreases  power  and  increases  throughput. 
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Network  Power  (dBW):  p  =- 1 16  dB 


Figure  30.  Power-rate  image  for  a  two-link  wireless  network  meeting  QoS  with  the  quadrant  marking  a 
Pareto  point. 
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5.  POWER-RATE  TRADE-OFFS 


Section  3  describes  the  geometry  of  the  two-link  mobile  wireless  network  operating  in  urban  scenar¬ 
ios.  This  section  describes  a  wireless  power  control  system  and  consequent  model  of  network  performance. 
Power  and  throughput  trade-offs  are  basic  performance  objectives  of  wireless  power  control  systems  oper¬ 
ating  in  power-limited,  ground-level  wireless  networks. 

Figure  3 1  is  a  schematic  of  a  generic  wireless  power  control  system  for  the  smallest  and  simplest 
two-link  mobile  wireless  network  of  Section  3.  Each  controller  commands  their  respective  transmitter  so 
that  the  received  SINR  “tracks”  the  specified  SINR.  These  controllers  are  in  conflict  because  the  trans¬ 
mit  power  on  one  link  is  the  interference  power  on  the  other  link.  This  model  of  a  wireless  power  control 
system  is  the  common  operational  framework  for  the  forthcoming  distributed  wireless  power  control  algo¬ 
rithms.  This  common  model  permits  comparisons  without  getting  entangled  in  system-specific  implemen¬ 
tations. 


PEJS) 


Figure  31 .  Two-link  schematic  of  a  generic  distributed  wireless  power  control  system. 

The  uplink  delays  and  downlink  delays  model  the  signal  processing  chain  to  implement  this  distributed 
wireless  power  control.  For  simplicity,  all  links  use  the  same  uplink  delays  and  the  same  downlink  de¬ 
lays.  The  downlink  delay  A tjj  models  the  processing  delays  mapping  the  analog  received  signal  into  the 
digital  controllers,  undertake  measurements  that  may  be  allowed  only  once  per  frame,  and  then  compute 
the  power  control  decision.  The  uplink  delay  At u  models  the  processing  delay  for  the  controllers  to  map 
the  power  control  decision  into  an  analog  power-control  message  that  is  broadcast  back  to  the  transmit¬ 
ters,  demodulating  this  message  at  the  transmitters,  and  then  implement  the  power  control  command  at 
the  transmitter’s  power  amplifier.  If  the  signals  are  narrowband,  the  wireless  power  flows  are  modeled  by 
Equation  (3)  and  repeated  here  for  convenience: 

P  i?(*)  =  G(t)pT(t)  +  Pjv(i). 

Figure  31  generalizes  the  wireless  power  flows  to  account  for  the  downlink  delay: 

PrW  =  G(t)pT(t  -  AtD)  +  p N(t). 
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Figure  3 1  also  shows  the  transmitters  setting  their  power  from  delayed  and  noise-corrupted  version  of  their 
receiver’s  power  command.  Although  the  wireless  system  describes  a  digital  messaging  scheme,  the  fol¬ 
lowing  model  eliminates  the  digital  processing  and  uses  the  power  control  vector  itself  to  set  the  transmit¬ 
ter’s  power: 

PtW  =  P c(*  -  Atu )  +  pB(i). 

The  wireless  power  control  models  in  this  report  lumps  both  delays  into  the  single  delay  r  and  omits  the 
error  vector.  Consequently,  the  transmitter’s  power  vector  is  modeled  as  delayed  version  of  the  control 
power  vector: 

PrO)  =  Pc(t~T)-  (8) 

Therefore,  power  control  in  the  wireless  network  is  the  narrowband  but  time-delayed  system: 

P *(*)  =  G(t)pc(t  -  t)  +  pN(t).  (9) 

The  update  commands  are  typically  issued  on  a  common  clock  on  each  link.  One  timing  scheme  resets  the 
transmitter’s  power  at  the  start  of  each  frame  in  the  received  signal.  Figure  32  illustrates  this  update-per- 
frame  and  the  timing — 1  millisecond  per  frame — used  in  this  report.  At  the  start  of  each  received  frame, 
receiver’s  controller  has  completed  the  power  measurements,  computed  the  power  control,  loaded  the 
power  control  into  a  message,  and  starts  the  broadcast  of  the  power-control  update  back  to  the  transmitter. 
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Figure  32.  Power  control  timing  with  the  update  issued  at  the  “downbeat”  of  each  frame. 


The  signal-to-interference  and  noise  ratio  (SINR)  bounds  the  throughput  of  the  network.  The  signal  power 
vector  is  measured  at  the  input  to  the  receiver’s  detector: 


P s(*)  =  Gs(t)pc(t  -  r);  Gs(t)  =  51 

The  interference  power  vector  is  measured  at  the  input  to  the  receiver’s 

P lit)  =  Gj(i)pcr(t  -  r);  Gj(t)  =  [  ° 

921 W 


0 

922  (t)_ 

detector: 

9t2  (*) 

0 


The  measured  SINR  vector  is 


7m(PcZ)  = 


P  s(t) 

Pjv(*)  +P/(i)' 


where  the  on  the  fraction  denotes  element-wise  division.  The  generalization  to  multiple  links  l  =  1, 
. . .  L  gives  the  standard  form  for  the  SINR  [79,  Equation  1]: 


_ 9u{t)pc,l{t  -  r) _ 

VN,l{t)  +  T,l^m9lm(t)PT,m{t  ~  t)  ' 


(10) 
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Figure  31  shows  the  controllers  compare  this  measured  SINR  7 M,i(f)  to  the  external  SINR  7a/,/ (i).  The 
controllers  command  their  transmitters  to  vary  their  broadcast  power  so  that  measured  SINR  “tracks”  the 
external  SINR.  The  tracking  error  and  the  adaptive  power  control  algorithms  arc  the  focus  of  this  report — 
not  the  problems  of  measuring  the  SINR. 

Measuring  SINR:  As  shown  in  Figure  31,  the  controllers  arc  limited  to  measurements  that  can  be  extracted 
from  the  received  signal  SRj(t)  as  seen  though  the  receiver’s  antenna  system.  Consequently,  measuring 
the  SINR  7 is  not  straightforward.  Most  radios  do  admit  measurement  of  the  total  received  power 
Piu(t).  Some  frames  allocate  slots  to  measure  the  noise  and  interference  p;.y,/  +  pij(t)  by  commanding 
the  transmitter  to  turn  off.  If  the  interference  dominates  the  noise,  7 admits  an  estimate  from  these 
measurements.  The  Kalman  filter  approach  recognizes  that  state-space  estimates  obtained  from  indirect 
measurements  arc  typically  higher  quality  than  downstream  approximations,  and  arc  obtained  with  less 
disruption  of  the  hardware  and  waveform  modifications.  For  example,  including  the  gain  gu  it)  as  paid  of 
the  state  allows  a  Kalman  filter  to  estimate  the  SINR  by  measuring  only  the  received  power  PRti(t). 

This  report  ranks  wireless  power  control  algorithms  assuming  each  controller  has  access  only  to  the 
SINR  7 The  assumption  allows  an  “apples-to-apples”  comparison  of  algorithm  performance  in  tac¬ 
tical  environments  without  the  confounding  effects  of  implementations.  Roughly  speaking,  optimizing 
network  performance  means  the  following: 


WPC-1  Minimize  network  power  py  (watts) 
WPC-2  Maximize  network  throughput  /A  (bps). 


Definition  7  (Network  Power)  The  network  power  is  defined  to  be  the  total  transmit  power: 

L 

Px(t)  ■■=  Y.PtM  [watts]. 

1=1 


Remarks :  (, i )  The  transmit  power  is  modeled  using  Equation  (8),  so  that  the  network  power  is  the  sum  of 
the  control  powers: 

L 

Pe(Pc;  t)  =  5>,i(f-r). 

i=i 

(;';')  This  network  power  model  does  not  include  the  power  required  for  messaging,  processing,  and  over¬ 
head. 

Definition  8  (Network  Throughput)  The  network  throughput  is  the  sum  of  the  link  capacities  integrated 
across  the  common  bandwidth 


L 

'■=  fB^2log2(l  +7M,l(Pc;*))- 
1=1 


Remarks:  (i)  The  network  is  narrowband,  (ii)  The  throughput  is  an  upper  bound  on  realizable  rates.  (Hi) 
Using  the  SINR  of  Equation  (10)  assumes  that  interference  is  independent  of  the  signal. 
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Definition  9  (Network  Performance)  The  network  performance  is  the  mapping  91  :  C7(M.  ®-+)  x  1  — > 
Ml  that  determines  the  instantaneous  trade-off  between  throughput  and  power. 


w(pc;*) 


Ps(Prfi£)l  watts 
,-Re(p  cl*).  Mbps 


(11) 


For  each  fixed  time  to,  the  network  performance  image  is  the  subset  of  the  power-rate  plane: 

9I(R+,t0)  :=  {9T(pc,t0)  :  Pc  G  M+}  ^  M+- 


Power  allocations  that  deliver  optimal  performance — maximum  throughput  and  minimum  power — lie  on 
the  upper  boundary  of  this  image.  At  these  Pareto  point,  where  no  change  in  power  can  simultaneously 
decrease  power  and  increase  throughput.  The  next  example  makes  this  geometry  explicit. 


Example  1  (Performance  Image)  Assume  the  gain  matrix  is  fixed: 


G 


TX-l  TX-2 


RX-l 

511 

512 

-72 

-75' 

RX-2 

.  521 

522 

-7! 

-80 

[dB], 


Set  the  noise  power  p^  =  —131  dBW.  Figure  33  shows  the  network  performance  image  for  the  transmit 
powers  sweeping —At)  <  pc.  i  ■  Pc. 2  <  0  dBW.  The  green  lines  correspond  to  constant  power  on  Link 
2.  The  blue  lines  correspond  to  constant  power  on  Link  1.  The  heavy  green  line  has  the  smallest  power 
allocated  to  Link  2.  The  heavy  blue  line  has  the  smallest  power  allocated  to  Link  1.  The  plot  shows  that  to 
maximize  throughput  and  minimize  power,  turn  Link  2  off  and  allocate  all  the  power  to  the  Link  1. 
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Figure  33.  Two-link  network  performance  image;  the  heavy  green  line  is  the  Pareto  front. 
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Figure  33  reports  that  the  optimal  network  performance  arises  by  allocating  all  network  power  to  Link 
1.  This  link  has  maximum  signal  gain.  However,  allowing  only  the  strongest  link  to  broadcast  voids  the 
meaning  of  “network”.  Consequently,  most  wireless  power  control  algorithms  aim  at  “fairer”  solutions  by 
enforcing  quality-of-service  (QoS)  constraints  [12],  [13].  One  common  QoS  requires  the  measured  SINR 
to  exceed  a  user- specified  external  SINR  on  all  links: 


7 x,i(t)  < 


_ gupc,i(t  ~  T) 

PN,l(t )  +  E  m^l  9lmPC,m(t 


t) 


7M,l(PCit ) 


for  l  =  1, . . . ,  L.  In  vector  notation. 


lx{t)< 


GsPcit  ~  T ) 

PnV)  +  GiPc(t  -  t) 


7m  (Pc!  t)i 


(12) 


where  the  “<”  denotes  the  vector  inequality  and  the  on  the  fraction  denotes  the  element-wise  division. 
That  is,  the  network  performance  is  to  be  optimized  under  the  constraint  that  each  link  is  guaranteed  a  min¬ 
imal  throughput. 


WPC-1  Maximize  network  throughput  /(v  (bps) 
WPC-2  Minimize  network  power  px  (watts) 
WPC-3  Guarantee  QoS. 


The  pitfall  of  this  approach  is  that  there  may  be  no  power  allocations  that  meet  the  QoS.  Consequently, 
all  approaches  to  the  power-rate-QoS  constraints  require  some  knowledge  of  the  feasible  control  vectors. 
Let 

VC(to)  ■=  {pc  G  :  7x(£o)<7m(Pc;*o)} 

denote  the  instantaneous  feasible  set  of  control  vectors  at  time  to-  The  network  performance  image  at  time 
to  is  the  subset  of  the  power-rate  plane: 

m(Vc,to)  ■■=  {7l(pc,t0)  :  Pc  G  Vc(to)}  C  K+. 

The  following  example  uses  the  same  gain  matrix  and  noise  level  as  in  the  proceeding  example.  The  QoS 
severely  constrain  the  resulting  network  performance  image. 

Example  2  (Performance  Image:  QoS)  Assume  the  noise  is  psj  =  —131  dBW  and  the  gain  matrix  is 


511 

<712 

'-72 

-75' 

521 

522 

-7! 

-80 

Figure  34  shows  the  resulting  network  performance  image  with  QoS:  7 m>  —  5  dB  or  equivalently,  each 
throughput  exceeds  0.395  Mbps.  Section  7  details  the  explicit  parameterization  of  this  power-rate  image. 

Figure  34  shows  the  power-rate  trade-offs  that  are  possible  under  the  QoS  constraint.  This  image  occu¬ 
pies  the  lower  edge  of  the  non-QoS  performance  image  in  Figure  33.  The  QoS  performance  image  shows 
that  massive  increases  in  network  power  have  little  effect  on  the  network  throughput.  Both  figures  point  to 
the  ideal  wireless  power  control  algorithm.  Such  an  algorithm  would  compute 
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Figure  34.  Two-link  network  performance  image  with  QoS.  The  upper  green  line  is  the  Pareto  front. 


•  the  Pareto  front  of  Vl(Vc,  to), 

•  the  control  vectors  that  parameterize  these  Pareto  fronts. 

This  distributed  wireless  power  control  is  computationally  challenging — the  multiobjective  optimization 
must  be  implemented  over  a  distributed  network  governed  by  the  delayed  control  system  of  Equation  (9). 
The  power-rate  trade-off  is  avoided  by  most  wireless  power  control  algorithms  by  constraining  one  of  the 
objectives.  Table  10  lists  selected  variations. 

Table  1 0.  Selected  single-objective  power  control  formulations. 


Description 

Objective 

Constraint:  0<p 

Power-constrained 
Rate-constrained 
Power-bounded 
Rate-power  Bounded 

max{f?s(p)} 

min{lTp} 

max{f?E(p)} 

min{lTp} 

Power:  =  1 1  p 

QoS:  Rb,x  <  Rb( P) 

Power:  p<pMAX 

QoS:  RbiX  <  Rb( p);  Power:  p<pMAX 

Nevertheless,  the  first  challenge  for  any  wireless  power  control  algorithm — distributed  or  centralized — 
is  to  determine  the  set  of  feasible  vectors.  Indeed,  the  QoS  constraint  may  be  so  stringent  that  no  power 
control  can  deliver  the  QoS.  Section  7  addresses  the  existence  of  feasible  power  control  Vc  0  for  the 
static  channel.  The  next  section  reviews  matrix  inequalities  to  facilitate  this  discussion. 


40 


6.  INTERLUDE:  MATRIX  INEQUALITIES 


Matrix  and  vector  inequalities  are  basic  tools  for  network  optimization.  There  arc  several  notations  for 
these  inequalities.  The  “<”  was  selected  for  similarity  to  the  “pointwise”  operations  of  MATLAB®. 

Definition  10  [46,  Chapter  15 ],  [47,  Chapter  16]  A  matrix  A  G  MMxjV  is  said  to  be  element-wise  non¬ 
negative  if  no  element  of  A  is  negative:  Am  n  >  0  or  A  >  0 . 

Just  as  scalar  inequalities  arc  preserved  under  multiplication  by  non-negative  scalars,  vector  inequali¬ 
ties  arc  preserved  under  multiplication  by  element-wise  non-negative  matrices. 

Lemma  3  [46,  page  531  ]  Let  A  £  WMxN  be  element-wise  non-negative:  A>0.  Ifx,  y  £  R;V  and  x<y, 
then  Ax<Ay. 

The  wireless  optimization  problems  reduce  to  constrained  optimization  problems  with  inequality  con¬ 
straints,  including  the  constraint  that  the  power  vector  p  is  non-negative: 

Definition  11  The  closed,  convex  orthant  of  M L  is  denoted  as 

:=  {p  £  Rl  :  p>0}  . 

Several  wireless  optimization  problems  reduce  to  linear  matrix  inequalities.  The  following  result  pa¬ 
rameterizes  a  matrix  inequality  when  the  matrix  is  invertible. 

Corollary  1  Suppose  A  £  MLxi  is  invertible.  Let  b  £  F®.  Then 

{x  £  Rn  :  Ax>b}  =  A_1b  +  A_1M+. 

Proof:  (  C):  The  inequality  Ax>b  is  equivalent  to  the  equality  Ax  =  b  +  y  for  some  y>0.  Apply  the 
matrix  inverse  A”1  to  the  equality  to  get  x  =  ,4 “  1  b  +  ,4“  'y.  (D):  Conversely  if  x  =  A~  b  +  ,4'  'y  for 
y>0.  Ax  =  b  +  y  or,  equivalently,  Ax>b.  Ill 

In  addition  to  the  element-wise  inequality,  there  is  the  element-wise  vector  product 

x  ©  y  :=  [xmym] 

and  the  element-wise  fraction 

x _ xm 

y  .  Vm.  _ 

Several  wireless  optimization  problems  require  inverting  a  matrix  of  the  form  I  —  F  where  F> 0.  If 
F  is  a  contraction,  the  geometric  series  computes  the  inverse  and  reveals  that  the  inverse  is  a  element-wise 
non-negative  matrix  because  each  term  Fk  is  an  a  element-wise  non-negative  matrix: 

OO 

{I-F)-1  =  ^Fk  >  0. 

k= 0 

The  key  assumption  is  that  F  is  a  contraction  relative  to  some  matrix  norm.  The  spectral  radius  determines 
this  matrix  norm. 
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Definition  12  (Spectral  Radius)  If  Ai,  . . . ,  A  l  are  the  eigenvalues  of  matrix  A  £  CLxL,  the  maximum 
extent  of  the  eigenvalues 

p(A)  :=  max{|A;|} 

is  called  the  spectral  radius  of  A. 

Theorem  2  [71,  page  384 ]  For  every  matrix  A  and  every  e  >  0,  there  exists  a  vector  norm  ||  o  | v  such  that 
the  induced  matrix  norm  is  e-bounded  by  the  spectral  radius: 

miU  :=  sup{Px||,  :  ||x||„  <  1}  <  p(A)  +  e. 

If  every  eigenvalue  A  of  A  with  |  A|  =  p(A)  has  only  linear  elementary  divisors,  there  is  a  vector  norm 
||  o  ||„  such  that  the  induced  matrix  norm  satisfies 

\\A\\v  =  p(A). 


Example  3  Ifb,  c  >  0,  the  spectral  radius  of 


is  the  geometric  mean  of  b  and  c:  p{F)  =  s/lx:.  If  po  >  0,  all  matrices  F  that  satisfy  p(F)  <  po  are 
parameterized  by  the  half-plane 

log(6)  +  log(c)  <  21og(p0)- 

If  F  has  a  spectral  radius  strictly  less  then  one:  p(F)  <  1,  there  is  a  matrix  norm  such  that  \\F\\  <  1 
and  ( I  —  F )  is  invertible  using  its  geometric  series  or,  equivalently,  the  power  method.  To  solve  (I  —  F) p  = 
y,  observe  the  part ial  sums 

Pa-  =  E  F y 

k= 0 

converge  to  the  solution  and  satisfy  the  recursion 

pk  =  y  +  fPk~i- 

The  parti al  sums  converge  according  to  the  tail  of  the  geometric  series  are  a  rate  determined  by  the  spectral 
radius. 


Lemma  4  (Partial  Sum)  Let  F  €  WLxL  with  \\F\\  <  1.  The  absolute  error  between  the  matrix  inverse 
(I  —  F  )  1  and  the  associated  partial  sums  is  bounded  as 


AbsErr(E;  K ) 


(I-F)-1 


K 


EF‘ 

fc=0 


F\\k+1 


n  =  0, 1, ... . 


Use  Theorem  2  to  select  a  matrix  norm 

ll-^ll  <  p{F)(l  +  e);  1  >  e  >  0. 

Lemma  4  bounds  the  absolute  error  as 

p(F)K+1(l  +  e)K+1 


AbsErr  (F;K)  < 


1  -  p(F) 
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7.  POWER-RATE  TRADE-OFFS  IN  STATIC  CHANNEL 

A  basic  problem  in  wireless  power  control  is  setting  the  minimum  SINR  requirements.  The  problem  is 
that  the  user-specified  SINR  yx  may  be  so  stringent  that  no  SINR  vector  7 m  exists  that  satisfies  7m  >7x 
or  that  the  set  feasible  control  vectors  Vc  is  empty. 

The  next  result  determines  when  feasible  control  vectors  exist — assuming  a  constant  channel,  a  con¬ 
stant  noise  power,  and  constant  'yx-  In  this  static  channel,  the  uplink  and  downlink  delays  arc  not  present 
and  the  set  of  feasible  control  functions  devolves  to  a  set  of  feasible  control  vectors. 

The  seminal  paper  in  wireless  power  control  was  written  in  1993  by  Foschini  &  Miljanic  [22].  This 
remarkable  paper  determined  that  the  set  of  feasible  control  vectors — if  nonempty — must  be  a  cone.  This 
cone  parameterizes  the  power-rate  image.  The  geometry  of  the  power-rate  image  determines  best  possible 
network  performance.  This  section  reviews  selected  results  of  the  Foschini  &  Miljanic  (FM)  paper  and 
presents  several  power-rate  images  showing  the  challenges  for  computing  optimal  performance. 


Theorem  3  (FM)  [22],  [75],  [82]  Let  G  £  RLxL  be  a  constant  gain  matrix  for  L  links.  Define  the  signal 
and  interference  matrices  as 


Gs  = 


9u  0  0 

0  522  0 

0  0  ••• 


0 

0 

0  9LL 


;  GI  =  G-GS. 


denote  the 


Let  pjy  £  denote  the  constant  and  positive  noise  vector:  Pjy>0.  kef  7 m  ■  ML  — » 
measured  SINR  vector 

(  \  _  GSPc 

lA1  Pc  Pn  +  giPc' 

Let  jx  denote  a  constant,  user-specified  SINR  constraint  vector.  Let  Vc  denote  the  set  of  feasible  control 
vectors: 

Vc  ■=  {Pc  €  RL  ■  1x<1m{Pc)}  ■ 

If  the  spectral  radius  of  the  feedback  matrix 

F  :=  Gg'TxGr,  Tx  :=  diag(7x) 


(13) 


is  strictly  less  than  unity:  p{F)  <  1,  the  vector 

Pmin  :=  (/  -  F^G^TxPn 

is  well-defined,  and  is  the  vertex  of  the  feasible  control  vectors: 

VC  =  Pmin  +  {!-  F^R1]  C  Rk. 
The  power-rate  mapping  91  :  Vc  C  — >  Ml  (Definition  9) 


N(Pc)  := 


Ps 

i?s 


watts 

Mbps 


(14) 


(15) 


(16) 


parameterizes  the  network ’s  power-rate  image  as 

m(Vc)  =  2R(pmin  +  (I-F)-1RI[) 
that  has  minimal  element:  9t(pmin)<91(T’c). 
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Proof:  Cross  multiply  to  get  the  linear  inequality: 


lx  < 


GsPc 

Pjv  +  GiPq 

lx  ©  (Pjv  +  GiPc)  <  GsVc 

^xPN  +  r YG/PC  <  G5PC;  r y  :=  diag(7x) 

TxPn  <  GsVc  —  FxGiPq 
Gs^xPn  ©  {I  ~  Gs'TxGfoc 
G^TxPn  <  (I  ~  F)Po 


(17) 


where  “©”  denotes  the  Hadamard  or  element-wise  product.  By  hypothesis,  p(F)  <  1.  Theorem  2  applies 
so  that  I  —  F  is  invertible  via  the  geometric  series 


(I  -F)-1  =  J2Fn  °' 

n= 0 

Lemma  3  applies  to  parameterize  all  vectors  pr;  <G  that  satisfy  Inequality  (17)  as  the  pointed  cone: 
{pc  G  ML  :  G^TxPn  <(/-  F) pc}  =  Pmin  +  (I  -  F)-1^. 


The  geometric  expansion  also  reveals  that  the  inverse  is  element-wise  non-negative  by  Equation  (13). 
Therefore,  pmin>0  and  the  set  of  vectors  that  satisfy  Inequality  (17)  must  also  lie  in  the  non-negative  or- 
thant: 

Pmin  + 

That  is,  all  the  vectors  satisfing  Inequality  (17)  arc  feasible  power  vectors.  For  the  minimal  element,  any 
feasible  vector  6  Vc  must  have  total  power  exceeding  the  minimal  power  lTpmin  <  lTpc.  by  Equa¬ 
tion  (15).  Likewise,  feasibility  forces 


lM  (Pmin)  =  lX  <1m(Pc)- 

Thus,  Tt(pmin)  is  the  minimal  element  of  the  power-rate  image.  Ill 

Three  observations  drop  out  of  this  result.  The  first  observation  is  that  this  result  determines  when  the 
set  of  feasible  vectors  is  non-empty — Vc  is  not  empty  provided  the  external  SINR  cx  is  sufficiently  small 
to  force  the  spectral  radius  p(  F)  <  1.  The  second  observation  is  that  the  gain  matrix  G  sets  an  upper  limit 
on  the  user-defined  SINR  or,  equivalently,  the  requested  network  throughput. 


Example  4  (External  SINR  Constraints)  The  narrowband  channel  gain  matrix  is 


The  matrix 


G 


TX-l  TX-2 


RX-l 

911 

912 

'-72 

-75' 

RX-2 

921 

922 

-7! 

-80 

[dB], 


F  :=  GfTxGr 


0  9nlx,i9i2 
922  lx:292\  0 


Tx  ■=  diag(7x) 


satisfies  the  inequality  p(F)  <  1  if  and  only  if  the  requested  SINR’s  are  bounded  by  the  excess  SINR  yc  as 
follows: 


lX,llX,2  < 


911922 

912921 


=■  IE- 
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Figure  35  illustrates  how  the  excess  SINR  77?  constrains  the  external  SINR’s  to  the  following  lower  half- 
plane: 

lx, 1  +  7X,2  <  (-72  -  80)  -  (-71  -  75)  =  -6  =  lE  [dB]. 

In  this  example,  positive  SINR  (dB)  cannot  simultaneously  exist  on  both  links  because  the  excess  SINR 
is  negative:  77?  =  —6  dB.  Figure  35  also  shows  that  the  gain  matrix  also  constrains  the  operational 
throughput — simultaneous  throughput  in  excess  of  0.6  Mbps  cannot  be  delivered  on  both  links. 


Figure  35.  External  constraints:  SINR  (left  panel),  throughput  (right  panel). 

For  networks  with  more  than  two  links  ( L  >  2),  the  spectral  radius  bound  p(F)  <  1  yields  an  upper 
bound  on  the  QoS: 

^  9ll 

7A7-  iiGzir 

This  upper  bound  has  the  physical  meaning  that  the  requested  SINR  7 x,i  cannot  exceed  the  channel’s  SIR  . 

The  third  observation  that  drops  out  of  Theorem  3  is  that  the  set  of  feasible  vectors  Vc  is  parameter¬ 
ized  by  the  closed,  convex  cone  under  the  affine  mapping: 

'Pc  =  Pmi„  +  (7  -  F)-1^. 

The  image  of  this  convex  cone  under  the  (nearly)  convex  power-rate  mapping  71  is  power-rate  (PR)  image 
for  the  static  channel: 

71(pc)=ffi(pmin  +  (/-F)-1M^). 

The  Pareto  front  of  this  power-rate  image  Tl('Pc)  is  the  best  possible  trade-off  between  power  and 
throughput  under  the  external  SINR  constraint. 


Example  5  (Power- Rate  Image)  The  narrowband  channel  gain  matrix  is 


G 


TX-l  TX-2 


RX-l 

511 

912 

-72 

-75 

RX-2 

921 

922 

-71 

-80 

[dB] 


from  Example  4.  Select  a  feasible  control  vector  less  than  the  maximal  SINR: 


+0 

-10 


< 


+3 

-9 
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lx  = 


[dB]. 


With  the  1-MHz  noise  set  to  psj  =  —116  clBW  and  rounding,  the  network  performance  image  91  (Pc)  of 
the  feasible  power  vectors 

PC'  =  Pmin  +  {I  ~  F)  1PK'i  P K  e 

Figure  36  shows  the  grid 'for  pK  £  (Left  Panel)  mapping  into  cone  of  feasible  vectors:  \>K  i-)-  pc 

(Right  Panel).  Figure  37  reports  how  this  grid  of  feasible  vectors  maps  into  the  network  performance  im¬ 
age:  pc  ^  «n(pc). 


Figure  36.  Pre-image  of  the  cone  of  feasible  power  vectors:  pK  e  M+  (left  panel).  Cone  of  feasible  power 
vectors:  pc  =  pmin  +  (J  -  F)~1pK ;  The  grid  is  distorted  by  the  dB  scaling  (right  panel). 


300-MHz;/  =  1  MHz  ;  0-dB  Coding  Gain 

1.9 


Figure  37.  Network  performance  image;  urban  noise  at  300  MHz  is  pn  =  -116  dBW;  the  green  dot  is 
91(Pmin)-  The  heavy  green  line  is  the  Pareto  front. 


46 


Mapping  the  cone  of  feasible  power  vectors  Vc  — >  Tt('Pc)  generates  the  network  image  in  Figure  37. 
This  image  shows  the  Pareto  front  is  a  function  of  the  total  network  power  for  this  two-link  network.  The 
next  result  verities  this  observation  is  true  for  networks  with  L  links  (Appendix  A  contains  the  proof). 


Theorem  4  (Pareto  Front)  Let  G  E  MixL  be  a  constant  gain  matrix  for  L  links.  Let  Gs  =  diag(G)  and 
Gi  =  G  —  Gs  denote  the  signal  and  interference  matrices,  respectively.  Let  p  v  E  denote  the  constant 
and  stricly  positive  noise  vector.  Let  7 m  ■  RL  — >  RL  denote  the  measured  SINR  vector 


7m(Pc)  =  ' 


GsPc 

Pjv  +  GiPc 


Let  yx  denote  a  constant,  user-specified  SINR  constraint  vector  with  Tx  =  diag(7x)  sufficiently  large  to 
force  the  feedback  matrix  F  =  Gf  ]T  xGj  to  have  the  spectral  radius  strictly  less  than  unity.  The  feasible 
control  vectors  are  the  positive  cone 


pc=pmm+(/-^r1Ri+cti 

with  vertex  pmin  =  (I  —  F)~1Gf1Tx  Pjv-  The  Pareto  front  of  the  network’s  power-rate  image  9T(' Pc)  is 
determined  by  maximizing  the  throughput  for  each  slice  of  the  total  network  power: 

ps  1  ^  max  {-Rs(Pc)  :  Pc  G  VcWY,  =  lTPc}  • 


The  next  example  shows  the  effect  of  increasing  noise  on  the  power  rate  image.  Roughly  speaking, 
increasing  noise  simply  shifts  the  image — the  Pareto  front  does  not  change. 


Example  6  (Power-Rate  Image  &  Noise)  The  narrowband  channel  gain  matrix  is 


G 


TX-l  TX-2 


RX-l 

5ii 

512 

-72 

-75 

RX-2 

521 

522 

-7! 

-80 

[dB] 


from  Example  4.  Select  a  feasible  control  vector  less  than  the  maximal  SINR: 


lx 


+0 

-10 


[dB]. 


Figure  38  shows  the  network  performance  image  when  the  noise  increases  to  px  =  —80  dBW. 


Example  6  also  shows  that  increasing  the  noise  power  only  shifts  the  performance  image  in  power.  The 
reason  is  that  measured  SINR  constraint  scales  with  the  noise: 


7m(Pjv!  p) 


Gsp 

Pjv  +  Gip 


If  p  G  meets  the  SINR  constraint:  7x<7a/(P.y:  p)  then  ap  meets  the  SINR  constraint  under  the  scaled 
noise: 


lx<  =  7jv/(apx;«p);  («  >  0). 

The  throughput  will  not  change  under  the  joint  scaling  of  the  power  and  noise: 

L 

Rb(  Pjv;p)  =  Ib  X/  log2  ( 1  +  7jw  (p  :V ;  p ) )  =  Rb(  pjv;p)- 


Z=1 
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Vt 

~ 

X 


f— 


0  2  4  6  8  10  12  14 

Network  Power  (dBW):  PN=- 80  dB 


300-MHz;/g=  1  MHz  ;  0-dB  Coding  Gain 


Figure  38.  Network  performance  image;  noise  power  pN  =  -80  dB;  the  green  dot  is  7t(prnin).  The  Pareto 
front  is  the  upper  green  line. 


Consequently,  the  network  performance  image  will  only  shift  in  power 


^(7A';aPiv) 


+  ^(7a;Pjv) 


with  the  noise  scaling — provided  the  power  is  plotted  in  dB.  However,  coding  gain  does  change  the  power- 
rate  image. 


Example  7  (Power- Rate  Image  &  Coding  Gain)  A  10-dB  coding  gain  is  modeled  by  reducing  the  off- 
diagonal  terms  by  a  factor  of  10.  Therefore,  the  narrowband  channel  gain  matrix  from  Example  4  is  modi¬ 
fied  as 


G 


TX-l  TX-2 


RX-l 

511 

512 

'-72 

-85' 

RX-2 

521 

522 

-8! 

-80 

[dB]. 


The  feasible  control  vector  is  less  than  the  maximal  S1NR  and  is  selected  to  be  the  same  as  the  preceding 
examples. 


lx 


+0 

-10 


< 


+13 

+1 


[dB]. 


With  the  1-MHz  noise  set  to  psj  =  —116  dBW,  the  network  performance  image 


K(Vc)  =  Vl(pmin  +  (I-FrlR2+) 


is  illustrated  in  Figure  39.  The  green  dot  is  71(pmin).  The  noise  power  in  each  link  is  pjz  =  —116  dBW. 
Figure  40  shows  the  corresponding  plot  when  the  noise  power  is  increased  to  p^  =  —80  dBW.  Compari¬ 
son  of  Figures  39  and  40  verifies  that  scaling  the  noise  simply  shifts  the  performance  image. 
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300-MHz;/fi=  1  MHz  ;  10-dB  Coding  Gain 


-38  -36  -34  -32 

Network  Power  (dBW):  /?  =- 1 16  dB 


Figure  39.  Network  performance  image  at  300  MHz  with  10-dB  coding  gain;  urban  noise  is  pN  =  -116  dB. 
The  Pareto  front  is  the  upper  green  line. 


300-MHz;/b=  1  MHz  ;  10-dB  Coding  Gain 


•2  0  2  4 

Network  Power  (dBW):  p  ~ 80  dB 


Figure  40.  Network  performance  image  at  300  MHz  with  10-dB  coding  gain;  noise  power  is  -80  dB.  The 
Pareto  front  is  the  upper  green  line. 
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The  final  example  extends  the  range  of  the  feasible  power  control  vectors  and  equalizes  the  SINR  re¬ 
quirement  to  —5  dB  on  both  links.  Figure  41  shows  that  the  network  mapping  “twists”  the  cone  on  the 
power  control  vector.  Standard  optimizers  seeking  the  Pareto  front  may  get  trapped  in  the  fold  and  con¬ 
verge  to  a  local  minimum. 


Gain  Matrix 
-72  -85 
-81  -80 
7X  (dB) 

-5 

-5 

Rx  (M-bps) 

0.396 

0.396 


300-MHz;  fg=  1  MHz  ;  10-dB  Coding  Gain 


■35  -30  -25  -20 

Network  Power  (dBW):  p  — 1 16  dB 


Figure  41.  Network  performance  image  at  300-MHz  with  10-dB  coding  gain;  noise  power  is  -116  dB;  the 
QoS  is  reduced  to  7 x  =  -5  dB  or  Rx  =  0.396  Mbps.  The  Pareto  front  is  the  upper  green  line. 
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8.  POWER-RATE  BASELINES  FOR  MOBILE  TWO-LINK  SCENARIOS 

Figure  42  is  a  top  view  of  the  suburban  Rosslyn,  VA,  showing  two  links: 

•  Link  1  is  fixed  with  both  TX- 1  and  RX-2  at  the  street  level 

•  Link  2  has  a  TX-2  traveling  along  the  .x-axis  track  while  RX-2  is  on  the  roof  of  an  80-meter  build¬ 
ing. 


Rosslyn,  Virginia 


Figure  42.  Top  view  of  Rosslyn,  VA,  showing  Link  1  (TX-1  to  RX-1)  and  Link  2  (TX-2  to  RX-2);  TX-2  travels 
along  the  x-axis  track. 


The  narrowband  gain  matrix  for  this  network  has  the  form 


G(x) 


RX-l 

RX-2 


TX-1  TX-2 

gn  gn(x) 

521  522  {x) 


[dB], 


where  the  first  column  is  constant.  Thus,  this  scenario  is  the  smallest  possible  mobile  network  to  bench¬ 
mark  the  wireless  power  control  algorithms. 

Section  8.1  shows  the  mobile  network  operating  along  the  track  at  300  MHz  using  1  watt  per  trans¬ 
mitter.  The  interference  from  Link  1  almost  shuts  down  Link  2.  The  mismatch  between  the  links  makes 
Quality-of-Service  challenging — the  minimal  throughput  over  the  25 -kHz  bandwidth  is  1 1  bps.  Sec¬ 
tion  8.2  shows  that  a  10-dB  coding  gain  mitigates  some  of  this  QoS  shortfall.  Both  transmitters  each 
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broadcast  1  watt  but  the  coding  gain  delivers  throughput  exceeding  25  kbps  over  63%  of  the  track.  The 
natural  question  asks:  what  is  the  payoff  of  using  wireless  power  control?  Section  8.3  determines  the  best 
possible  payoff  for  wireless  power  control  using  the  following  assumptions: 

•  Centralized  control 

•  Channel  matrix  is  perfectly  known 

•  Noise  power  is  perfectly  known 

•  No  delay. 

That  is,  this  network  is  equipped  with  an  all-knowing  “genie.”  This  genie  is  tasked  to  allocate  power  to 
maximize  throughput  while  maintaining  a  minimal  QoS  of  25  kbps.  Best  possible  power  allocation  boosts 
this  QoS  to  cover  74%  of  the  track.  Section  8.4  moves  the  frequency  region  to  1.9  GHz  and  increases  the 
bandwidth  to  1  MHz. 


8.1  MOBILE  TRACK  AT  300  MHZ:  FIXED  POWER 


Figure  43  plots  the  gain  matrix  generated  by  TX-2  moving  along  the  x  axis: 


G{x) 


TX-l 

TX-2 

511 

512  (z) 

'-72 

512  (®) 

521 

522  O)  . 

-71 

522  0)_ 

[dB], 


These  gains  exhibit  multiple  phenomena:  the  characteristic  lobing  when  only  a  few  strong  rays  interfere, 
the  sharp  jumps  when  a  street  “unblocks”  a  transmitter,  and  the  large-scale  trends  caused  by  distance. 


Rosslyn  at  300  MHz:  Moving  TX-2;  Gain  Matrix 


-115 


0  50  100  150  200 

x  (meters) 


Figure  43.  Gain  matrix  for  moving  TX-2  traveling  along  the  x  axis. 
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The  signal  gain  for  Link  1  falls  below  the  interference  for  10  <  x  <  75:  gn  =  —72  dB;  fj\2(x)  > 
—70.  If  equal  transmit  powers  are  employed  and  the  noise  is  relatively  small.  Link  1  is  interference -limited 
over  this  region  and  has  negative  SINR.  The  signal  gain  for  Link  2  is  comparable  to  interference  over  the 
whole  track:  <722  ~  (J2\  {■>')■  If  equal  transmit  powers  are  employed  and  the  noise  is  relatively  small.  Link  2 
will  have  a  relatively  flat  SINR  around  0  dB.  Figure  44  reports  the  SINR  produced  by  moving  TX-2  along 
the  x-axis  assuming: 

•  Unit  power:  pr  =  [1;  1]T  watts 

•  Low  noise:  pn  =  —133  dBW. 


Rosslyn  at  300  MHz:  Moving  TX-2 


-30 


0  50  100  150  200 

x  (meters):  p  =[1;  1]  Watts 

Figure  44.  SINR  for  moving  TX-2  traveling  along  the  x-axis. 

For  Link  1,  the  power  from  TX-1  to  RX-1  is  not  changing,  so  the  “S”  in  SINR  is  constant,  whereas  the 
power  from  TX-2  is  rolling  off  as  TX-2  travels  along  the  x-axis  (<712(3:)  |).  Consequently,  RX-1  receives 
less  interference  from  TX-2  and  shows  an  increase  in  the  measured  SINR  7m, t-  For  Link  2,  the  power 
from  TX-1  to  RX-2  is  not  changing,  so  the  “I”  in  SINR  is  constant,  whereas  the  power  from  TX-2  is  in¬ 
creasing  as  TX-2  travels  along  the  x-axis  (<722(3:)  f).  Consequently,  RX-2  receives  more  signal  from  TX-2 
and  also  shows  an  increase  in  the  measured  SINR  7m, 2- 


53 


Figure  45  reports  the  throughput  produced  by  moving  TX-2  along  the  x  axis  assuming: 

•  Unit  power:  pr  =  [1;  1]T  watts 

•  Low  noise:  pxt  =  —133  dBW 

•  Bandwidth:  fs  =  25  kHz. 


Rosslyn  at  300  MHz:  Moving  TX-2;  f  =  25  kHz 

B 


0  50  100  150  200 


T 

x  (meters):  p  =[1;  1]  Watts 

Figure  45.  Throughput  for  moving  TX-2  traveling  along  the  x  axis. 

The  throughput  is  totally  dominated  by  Link  1  because  the  SINR  of  Link  1  typically  exceeds  the 
SINR  of  Link  2  by  10  dB.  Giving  equal  power  to  both  links  “throws  away”  total  throughput  because  of 
the  SINR  disparity  between  the  links.  Therefore,  throughput  would  be  maximized  if  all  the  power  was 
allocated  to  Link  1  and  Link  2  was  shut  down.  If  a  minimal  quality-of-service  (QoS)  is  enforced,  through¬ 
put  is  reduced  because  the  throughput  in  Link  1  is  traded  to  boost  throughput  in  Link  2.  That  is,  the  mea¬ 
sured  SINR’s  must  exceed  the  external  SINR:  7m >7 v-  Indeed,  Figure  45  uses  an  implicit  external  SINR: 
7x  =  —  35  [dB]  set  by  Figure  44.  The  QoS  is  equivalent  to  throughput: 

Rb  =  fs  log2(l  +  7)  5=>  7  =  2Rb/fB  -  1. 

Therefore,  xx  =  —35  dB  implies  an  equivalent  QoS  Rx  =  11.4  [bps],  assuming  a  bandwidth  Jb  =  25 
kHz.  Coding  gain  can  mitigate  this  low  QoS. 
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8.2  MOBILE  TRACK  AT  300  MHZ:  CODING  GAIN 


This  section  applies  a  10-dB  coding  gain  to  the  gain  matrix: 


G(x) 


TX-l  TX-2 


RX-l 

911 

9cl9i2{x) 

1 

to 

9i2(x)  -  10 

RX-2 

.  9c1 921 

922  (x) 

-8! 

922  (x) 

[dB] 


to  get  a  network  with  practical  link  values.  Figure  46  reports  the  SINR  produced  by  moving  TX-2  along 
the  x  axis  assuming: 


•  Unit  power:  pr  =  [1;  1]T  watts 

•  Fligh  noise:  pn  =  —80  dBW 

•  Coding  gain:  gc  =  10  dB. 


Rosslyn  at  300  MHz:  Moving  TX-2;  Coding  Gain=10  dB 


SINR:  Link  2 


_ _ i _ I _ ! _ I _ I _ 

0  50  100  150  200 

x  (meters):  p  =[1;  1]  Watts 

Figure  46.  SINR  for  moving  TX-2  traveling  along  the  x  axis;  10-dB  coding  gain. 

The  first  observation  is  that  SINR  on  both  links  are  comparable  because  of  the  coding  gain.  The  sec¬ 
ond  observation  is  that  “flat”  areas  for  the  SINR  on  Link  1  are  caused  by  the  coding  gain  and  the  increased 
noise:  pn  =  —80  dBW.  The  noise  was  increased  because  of  the  current  limitations  on  the  Pareto-front 
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code.  Figure  43  showed  that  Link  1  was  interference  limited  for  10  <  x  <  75  (meters).  Exterior  to  this 
region,  the  interference  is  smaller  than  the  noise:  Over  this  region, 

9c1  9i2(x)pt,2  =  -10  -  100  +  0  <  -80  =  pn  [dB], 

In  this  low-interference  region,  the  SINR  on  Link  1  is  almost  constant  as  shown  in  Figure  46: 

7M,i(z)  = - 9itPt,i  - _  SRpT1  =  _72  +  80  +  0  =  8  [dB] 

Pn  +  9C  9i2[x)pt,2  Pn 

Figure  47  reports  the  throughput  determined  by  the  noise  and  coding  gain.  The  figure  suggests  that  the 
large  throughput  on  Link  1  can  be  traded  to  boost  the  smaller  throughput  on  Link  2  so  that  a  minimal  rate 
of  25  kbps  is  possible. 


Rosslyn  at  300  MHz:  Moving  TX-2;  f  =  25  kHz 

B 


0  50  100  150  200 


T 

x  (meters):  p  =[1;  1]  Watts 

Figure  47.  Throughput  for  moving  TX-2  traveling  along  the  x  axis;  Rosslyn  at  300  MHz;  1 0-dB  coding 
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8.3  MOBILE  TRACK  AT  300  MHZ:  BEST  PERFORMANCE 


Figure  48  reports  the  throughput  for  this  noise,  coding  gain,  and  the  explicit  QoS  constraint: 

Rb,x  =  25  [kbps]  =  0  [dB] . 

The  figure  reports  “Link”  increased  from  63%  to  74%.  The  percentages  arise  when  both  Link  1  and  Link  2 
exceed  the  explicit  QoS  constraint  over  the  track.  The  first  number  comes  from  Figure  46  where  the  power 
was  equally  divided  between  both  links.  This  implicit  design  delivered  25-kbps  links  over  63%  of  the 
track.  The  second  number  is  delivered  by  a  Pareto-front  optimizer  that  distributes  the  the  power  between 
the  two  links.  Specifically,  the  network  “genie”  seeks  to  maximize  throughput  by  allocating  total  power 
while  guaranteeing  QoS: 

{L  ' |  (  Positive  powers:  p>0 

i?s(p)  =  ^  Rb,i(p)  z  subject  to  <  Fixed  network  power:  Vt,  =  Yli=\Pi  • 

1=1  '  1  QoS:  Rb,i(p)>Rb,x 

This  is  the  power-constrained  optimization  of  Table  10.  Figure  48  is  produced  by  applying  this  Pareto- 
front  optimizer  to  each  gain  matrix  along  the  track.  This  best  possible  centralized  control  delivered  the  25- 
kbps  links  over  74%  of  the  track.  Little  difference  is  observed  between  this  optimized  throughput  and  the 
equal-power  design — except  the  parts  of  the  track  where  the  Pareto-front  optimizer  delivers  the  improved 
QoS. 


Rosslyn  at  300  MHz:  Moving  TX-2;  f  =  25  kHz 

B 


20 

0 


0  50  100  150  200 


x  (meters):  Total  Power=2  (Watts);  QoS>  25  kbps;  Link:  63— >74% 


Figure  48.  Throughput  for  moving  TX-2  traveling  along  the  x  axis;  10-dB  coding  gain;  QoS>25  kbps.  The 
“blanks”  in  the  plot  mark  where  QoS  failed. 
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8.4  MOBILE  TRACK  AT  1 .9  GHZ 


Figure  43  plots  the  gain  matrix  produced  by  moving  TX-2  along  the  x  axis  operating  at  1.9  GHz: 


G(x) 


-102  512  (x) 
-93  522  (x) 


[dB], 


In  comparison  to  the  propagation  at  300  MHz  as  seen  in  Figure  43,  the  channel  matrix  has  the  same  gen¬ 
eral  form:  RX-1  is  still  jammed  by  TX-2  at  the  beginning  of  the  track  and  then  is  blocked  by  buildings 
once  TX-2  moves  past  x  =  75  meters;  there  is  the  characteristic  “2 -ray”  lobing  on  both  522  (x)  and  512  (x) 
at  the  end  of  the  track. 


Rosslyn  at  1900  MHz:  Moving  TX-2;  Gain  Matrix 


PQ 

-cs 


-120 


100 

x  (meters) 


Figure  49.  Gain  matrix  for  moving  TX-2  traveling  along  the  x  axis;  Rosslyn  at  1 .9  GHz. 
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Figure  50  shows  the  SINR  on  both  link  in  a  high-noise  environment.  Suburban  noise  is  typically  p^  = 
—  149  dBW  at  1.9  GHz  over  the  fs  =  1.3-MHz  band.  These  simulations  use  pjy  =  —100  dBW  for  a 
high-noise  environment.  Link  1  shows  the  flat  regions  characteristic  of  high  noise:  the  noise  dominates  the 
interference  ( gi2PT,2  -C  pm)  so  that 


_  fj\  \Pt,\ 

PN  +912(x)pt,2 
~  9uPt,i 

PN 

=  gu+PT,l-PN  [dB] 
=  -102  +  0  +  100  [dB] 
=  -2  [dB], 


These  flat  regions  will  control  the  link  QoS. 


Rosslyn  at  1 .9  GHz:  Moving  TX-2;  Coding  Gain=10  dB 


SINR:  Link  1 


_,n .  . SINR:  Link  2 

_ i _ i _ i _ L _ i _ 

0  50  100  150  200 

T 

x  (meters):  p  =[1;  1]  Watts 

Figure  50.  SINR  for  moving  TX-2  traveling  along  the  x  axis;  Rosslyn  at  1 .9  GHz;  coding  gain  gc  =  10  dB; 
equal  power  pTtl  =  pt, 2  =  1  watt. 
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Figure  5 1  shows  the  throughput  for  this  high-noise  simulation.  The  network  throughput  is  dominated 
by  Link  2  because  Link  1  is  limited  by  the  high  noise:  i2&,i  (x)  <  1  Mbps  in  the  noise-limited  regions. 


Rosslyn  at  1.9  GHz:  Moving  TX-2;  f  =  1.3  MI Iz;  Coding  Gain=T0  dB 
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Figure  51 .  Equal-power  throughput  for  moving  TX-2  traveling  along  the  x  axis;  Rosslyn  at  1 .9  GHz;  coding 
gain  gc  =  10  dB;  equal  power  pT,i  =  pt, 2  =  1  watt. 


Figure  52  illustrates  best  possible  power  control  the  QoS  requirement — each  link  carries  at  least 
Rt„,x  =  1  Mbps.  This  best  possible  control  assumes  a  network  “genie”  maximizing  throughput  subject 
to  QoS  constraints.  At  each  point  on  the  track,  the  genie  solves 


max 


^2Rb,l(  P) 

z=i 


subject  to 


Positive  powers: 
Fixed  network  power: 
QoS: 


P>0 

PS  =  Ez= 1  PI 
Rb,i(p)>Rb,x 
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In  the  equal  power  case  of  Figure  51,  Link  1  never  meets  the  QoS:  Rbyi(x)  <  Rb,x  so  that  QoS  is 
0%  along  the  track.  However,  Figure  5 1  suggests  that  some  of  the  throughput  on  Link  2  can  be  traded  to 
boost  the  throughput  on  Link  1  to  meet  the  QoS.  Figure  52  reports  that  this  trading  is  possible — QoS  is 
delivered  over  54%  of  the  track.  The  plot  shows  where  the  QoS  is  delivered  by  plotting  both  links  and  their 
total.  The  blank  regions  mark  where  the  QoS  requirements  cannot  be  met. 


Rosslyn  at  1.9  GHz:  Moving  TX-2;  f .=1.3  MHz;  Coding  Gain=10  dB 
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J_ I_ I_ I  I 

0  50  100  150  200 

x  (meters):  Total  Power=2  (Watts);  QoS>  1  Mbps;  Link:  0— »54% 


Figure  52.  Throughput  for  moving  TX-2  traveling  along  the  x  axis;  Rosslyn  at  1 .9  GHz;  coding  gain 
gc  =  10  dB;  QoS  >  Rb,x  =  1  Mbps.  The  “blanks”  in  the  plot  mark  where  QoS  failed. 
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9.  THE  FOSCHINI  &  MILJANIC  (FM)  ALGORITHM 


In  1993,  Foschini  &  Miljanic  produced  their  seminal  paper  on  distributed  wireless  power  control  [22]: 
A  Simple  Distributed  Autonomous  Power  Control  Algorithm  and  Its  Convergence.  Their  excellent  solution 
characterized  the  existence  of  solutions  (feasibility),  introduced  the  concept  of  the  all-knowing  network 
“genie,”  computed  the  global  solution  via  their  network  “genie”  inverting  a  matrix,  delivered  a  distributed 
solution  based  on  a  simple  feedback  ODE  whose  discrete  version  is  equivalent  matrix  inversion  via  the 
power  method,  and  finally,  analyzed  convergence  and  stepsize  control.  Assume  th e.  feedback  matrix  of 
Equation  (13) 

F  =  G^TxGi;  Tx  '.=  diag(7x) 

has  spectral  radius  strictly  less  than  unity:  p(F)  <  1.  The  minimum  power  allocation  pmin>0  that  meets 
the  QoS  7M(pmin)  =  lx  is  [22,  Equation  12]: 

Pmin  :=  ~  FT'Gg'TxVN- 


The  power  method  for  matrix  inversion  produces  the  sequence  of  power  control  vectors  \i(  it)  at  update 
times  tn  [22,  Equation  14]: 


Pmin  F-  pc(tn+l) 


Gs^xPN  +  Fpitn) 
G/TxPx  +  G/TxGjpdtn) 
G^Tx  {Pjv  +  Gjpc(tn)} 

G^Tx  {Pjv  +  P i{tn)}  ■ 


The  power  method  is  distributed  along  each  link  l  =  1, . . .  L  because  [22,  Equation  2]: 

Pmin,/  F-  PC,l{tn+ 1)  =  9^1X1  {PN,l  +  Pl,l(tn)} 


=  PC,l(tn)lX,l 


PN,l  +Pl,l{tn) 
9llPC,l(tn ) 


=  PC, l  (tn) 


IX, l 
lM,l(tr. 


That  is,  each  receiver  only  has  to  measure  its  SINR  to  implement  the  update.  If  the  wireless  network  has  a 
system  delay,  the  receiver  measures  the  delayed  SINR: 

,  =  _ 9ll(tn)PC,l(tn  -  T) _ 

lM)l  U  PN,l{tn)  +  9lm(tn)pC,m(tn  ~  t)  ' 

The  following  algorithm  includes  a  stepsize  (3  to  control  stability  using  these  time-delayed  updates. 


Algorithm  1  (FM)  [22]:  [27],  [79]  Select  the  step  size  0  <  /3  <  1. 

FM-1  The  network  is  synchronized  and  all  updates  occur  at  the  star t  of  each  frame:  tn  :=  ntp. 

FM-2  Each  receiver  measures  its  SINR  at  the  start  of  each  frame: 

(t  S  =  _ 9ll{tn)PC,l{tn  -  t) _ 

lM'1  n  PN,l(tn)  +  9lm  (fn)PC,m  ifn  ~  t)  ' 
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FM-3  The  power  is  updated  for  the  frame  [tn,  tn+ 1)  as 

PcAtn+ 1)  =  PC,l(fn)  {  7A/  ,  | 
l  TM,l  (fn  )  J 

FM-4  The  control  power  is  limited 

PC,lifn+ 1)  =  [Pmin  5;  PC,l(^n+ 1)  tS  Pmax] • 


9.1  FM  TRACK  PERFORMANCE  AT  1.9  GHZ 

Table  1 1  lists  the  RF  parameters  of  the  gain  matrix,  the  SINR,  and  the  updating  of  the  wireless  control 
algorithms.  These  parameters  are  outside  the  control  of  the  FM  algorithm. 

Table  11.  RF  and  system  parameters. 


Variable 

Value 

Description 

fc 

1.9  GHz 

Carrier  frequency 

fs 

1.3  MHz 

Bandwidth 

Pn 

-149  dBW 

RF  noise 

T 

0.9  ms 

System  delay 

9c 

10  dB 

Coding  gain 

tF 

1  ms 

Frame  duration 

Table  12  lists  the  inputs  to  the  FM  algorithm.  The  QoS  is  1  Mbps  on  both  links;  the  frame  duration 
t-F  =  1  ms  sets  the  power  control  update  rate. 

Table  12.  FM  parameters. 


Variable 

Value 

Description 

Rb.x 

Pmin 

P  max 

Pini 

P 

1  Mbps 

0  W 

1  W 

1 0  /iW 
0.9 

Throughput 
Minimum  power 
Maximum  power 
Initial  power 
Stepsize 

Figure  53  reports  the  power  and  the  rate  for  the  network  during  the  time  TX-2  travels  along  the  x-axis 
track.  The  upper  panel  compares  the  power  on  each  link  computed  by  the  distributed  FM  algorithm.  The 
lower  panel  compares  throughput  on  each  link.  Although  the  coding  gain  is  gc  =  10  dB,  throughput  is  set 
using  the  capacity:  Rif  t)  =  fn  x  log2(l  +  7m(0)-  On  the  average,  the  FM  algorithm  delivers  rates  at 
Rb,x  =  1  Mbps. 

Figure  54  assesses  the  average  performance  of  the  FM  algorithm  as  the  mobile  network  travels  along 
the  x  axis.  Also  plotted  is  an  upper  bound  on  performance  given  by  the  Exact  FM  algorithm  (description 
forthcoming).  Neither  algorithm  meets  the  QoS  constraints  at  every  point  of  the  x-axis  track.  Therefore, 
the  next  two  paragraphs  detail  how  the  average  performances  account  for  this  loss  of  QoS. 

FM  Algorithm:  let  Qfm  index  those  frames  meeting  the  QoS: 

Qfm  :=  {n  £  [1  :  Np]  -  1  :  7Af(Pcl  tn)> 7.\7  tn  =  nip}, 
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Rosslyn  at  1.9  GHz:  TX-2=4  m/s;  Py=- 149  dBW 
FM  Algorithm:  Coding  Gain=10  dB;  /?=0.9 


Figure  53.  FM  performance  for  mobile  TX-2  traveling  along  the  x  axis;  Rb>x  =  1  Mbps;  delay  r  =  0.9  ms; 
frame  duration  t,F  =  1  ms;  fc  =  1.9  GHz;  Coding  gain  gc  =  10  dB.  Upper  Panel:  Link  power  along  the 
track;  Lower  Panel:  Link  rates  along  the  track. 


where  Np  is  the  number  of  frames  along  the  track.  Averaging  the  rates  over  the  entire  track  mixes  the  un¬ 
successful  QoS  transmission  with  the  transmission  that  met  the  QoS.  This  mixing  is  problematic  when  one 
link  dominates.  Indeed,  Section  8.4  reports  Link  1  always  delivers  high  rates  although  Link  2  was  almost 
shut  down — total  throughput  is  high  while  the  QoS  is  never  satisfied.  Therefore,  the  average  network  rate 
is  defined  as  average  rate  only  when  the  network  delivers  QoS: 


Rfm 


1 

Wp 


E  iT 

neQFM 


Rb,l(tn) 

Rb,  2  (J'n) 


The  FM  algorithm  always  transmits  minimum  power  regardless  of  the  QoS.  Therefore,  network  power  is 
the  average  power: 

1  nf- i 

PFM  v  X!  lTP  din)- 

F  71=0 


Exact  FM  Algorithm:  This  algorithm  operates  with  zero  system  delay,  global  network  knowledge  (i.e.,  the 
gain  matrix),  and  solves  for  minimum  power  pmin  meeting  the  QoS  constraints.  Not  every  position  on  the 
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Rosslyn  at  1.9  GHz:  TX-2=4  m/s;  p  =- 149  dBW 
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Figure  54.  FM  performance  averaged  over  x-axis  track;  QoS  Rbtx  =  1  Mbps;  delay  r  =  0.9  ms;  frame  rate 
tF  =  1  ms;  fc  =  1.9  GHz;  coding  gain  gc  =  10  dB. 


.x-axis  track  admits  a  solution.  Let  Qxfm  index  those  frames  meeting  the  QoS: 

Qxfm  :=  {n  G  [1  :  Np]  —  1  :  7M(Pminfa))>7x;  tn  =  ntp}. 
Define  the  average  network  rate  as  the  average  rate  only  when  the  network  delivers  the  QoS: 


1 


RxFM=m  ^  ^ 


Rb,x 

_Rb,  x 


=  2 Rb,x  x 


neSxFM 

The  average  power  used  by  the  exact  FM  algorithm  is  defined  as 

1 


Qxfm 
N,  ' 


PXFM  = 


I  Qxfm  | 


lTPmin(*n). 


neQxFM 

Indeed,  the  power  is  not  defined  when  the  QoS  constraint  cannot  be  satisfied. 


Discussion :  Figure  54  reports  that  the  FM  algorithm  delivers  approximately  70%  of  the  best  possible  av¬ 
erage  network  rate  using  only  slightly  more  power.  The  FM  delivers  this  performance  against  the  system 
delay  of  r  =  0.9  ms  and  no  global  knowledge  of  the  network — each  link  uses  only  information  measured 
at  its  receiver.  The  next  section  assesses  the  FM  algorithm  against  larger  system  delays. 
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9.2  FM  PERFORMANCE  AGAINST  DELAY 


Figure  55  compares  the  average  power-rate  performance  for  the  Exact  FM  algorithm  to  the  distributed 
FM  algorithm.  As  in  the  Section  9.1,  the  exact  FM  algorithm  operates  with  zero  system  delay,  global  net¬ 
work  knowledge,  and  solves  for  minimum  power  pmin  that  meets  the  QoS  constraints.  Therefore,  its  per¬ 
formance  has  not  shifted  from  Figure  54.  The  distributed  FM  algorithm  “hangs”  tight  until  the  delay  ex¬ 
ceeds  8  ms.  At  9  ms  delay,  the  power  jumps  by  15  dB.  Increasing  delays  in  excess  of  9  ms  force  the  power 
to  increase  while  the  throughput  rolls  off.  Thus,  the  FM  algorithm  has  a  delay  threshold — exceeding  this 
threshold  incurs  system  collapse. 
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Figure  55.  Average  FM  performance  along  i-axis  track  as  a  function  of  delay;  RbtX  =  1  Mbps;  frame  rate 
tF  =  1  ms;  fc  =  1-9  GHz;  coding  gain  gc  =  10  dB. 
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10.  THE  MODPC  ALGORITHM 


An  excellent  generalization  of  the  FM  algorithm  is  the  Multiple  Objective  Distributed  Power  Control 
Algorithm  (MODPC)  of  Elumsrati,  Jantti,  and  Koivo  [20].  Their  generalization  starts  with  the  observa¬ 
tion  that  the  FM  algorithm  seeks  to  minimize  power  subject  to  the  QoS  constraint.  The  MODPC  algorithm 
[20]: 


“  . . .  treats  the  QoS  requirement  an  another  objective  for  the  power  control  and  [produces]  a 
fully  distributed  method  for  solving  the  multiobjective  power  optimization  problem.” 


Their  derivation  starts  by  weighting  the  excess  power  and  the  QoS  gap  [20,  Equation  3]: 


ei(t)  =  \p\pc,l(t )  -  Pmin, 1 1  +  Xq\ix,1  ~  lM,l(t)\ 


on  each  link  l  =  1,  . . . ,  L  where  the  power  weight  Xp  >  0  and  the  QoS  weight  Xp  >  0  have  unit  sum: 
Xj>  +  Xq  =  1.  The  global  cost  function  sums  these  errors  over  the  entire  network  and  a  time  window  [20, 
Equation  5] : 

L  M 

J{PC^n)  =  ^  ^  ^  ® p  e-iitn—m)  • 

/=1  m= 0 

Here  the  time  window  looks  back  M  time  slots  and  weights  the  past  with  the  “forgetting”  factor  up  >  0. 
A  power  control  is  sought  that  minimizes  this  global  cost  function.  Under  the  1-tap  AR  assumption  of  the 
power  evolution  [20,  Equation  6],  the  MODPC  power  control  is  [20,  Equation  7] 


PC,l(tn+ 1) 


PC,l(tn ) 


ApPmin  +  Xq'Jx 
XpPC,l{tn )  +  XQ^M,l(tn) 


This  solution  contains  the  FM  algorithm  for  Xp  =  0: 


PC,litn+ 1) 


Pc, i  (tn) 


lx 

lM,l{tn) 


Our  interest  lies  at  the  other  extreme  when  Ap  «  1.  Elumsrati,  Jantti,  and  Koivo  [20]  analyzed  conver¬ 
gence,  and  compared  their  MODPC  algorithm  against  the  DPC  algorithm  [23],  FM  algorithm  [22],  and 
a  second-order  power  control  (SOPC)  algorithm  [35],  The  remainder  of  this  section  ranks  the  following 
MODPC  algorithm  against  the  exact  FM  algorithm. 


Algorithm  2  (MODPC)  Select  the  step  size  0  <  /3  <  1  and  non-negative  weights  X p  and  Xq  such  that 
Xp  +  Xq  =  1. 


MODPC-1  The  network  is  synchronized  and  all  updates  occur  at  the  start  of  each  frame:  tn  :=  ntp. 
MODPC-2  Each  receiver  measures  its  SINR  at  the  start  of  each  frame: 

(  ,  =  _ 9ll{tn)PC,l{tn  -  T) _ 

lM)l  n  PN,l  (tn)  +  Em/I  9lm  ( tn)PC,m  ( tn  ~  t)  ' 

MODPC-3  The  power  is  updated  for  the  frame  [tn ,  tn+  \ )  as  [20,  Equation  7] 

/,  \  / ,  \  /  Xppmjn  +  Xqjx  1 

PC,l(tn+ 1)  =  PC,l(tn)  <  T - TTTTA - 7TX  I 

{  Xppc,i(tn)  +  Xqim,i (tn )  J 

MODPC-4  The  control  power  is  limited 

PC,l(tn+ 1)  =  [Pmin  A  PC,l(tn+ 1)  A  Pmax]  • 


1-/3 
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10.1  MODPC  TRACK  PERFORMANCE  AT  1.9  GHZ 


Table  1 1  lists  the  RF  and  system  parameters  common  to  all  the  algorithms.  These  parameters  are  out¬ 
side  the  control  of  the  MODPC  algorithm.  Table  13  lists  the  inputs  to  the  MODPC  algorithm.  The  QoS 
is  1  Mbps  on  both  links;  the  frame  duration  tp  =  1  ms  sets  the  power  control  update  rate.  The  table  sets 
MOPDC  algorithm  far  from  the  FM  algorithm  by  selecting  the  weights  A p  =  0.99  and  A q  =  0.01,  This 
weighting  minimizes  power  at  the  expense  of  QoS. 

Table  13.  MODPC  parameters. 


Variable 

Value 

Description 

Rb,X 

1  Mbps 

Throughput 

Pmin 

0  W 

Minimum  power 

Pmax 

1  W 

Maximum  power 

Pini 

1 0  /jW 

Initial  power 

P 

0.9 

Stepsize 

A  p 

0.99 

Power  weight 

XQ 

0.01 

QoS  weight 

Figure  56  reports  the  power  and  the  rate  for  the  network  during  the  time  TX-2  travels  along  the  x-axis 
track.  The  upper  panel  compares  the  power  on  each  link  computed  by  the  MODPC  algorithm.  The  lower 
panel  compares  throughput  on  each  link.  Although  the  coding  gain  is  gc  =  10  dB,  throughput  is  set  using 
the  capacity:  Rb(t)  =  fn  x  log2(l  +  7 M{t)).  On  the  average,  the  MODPC  algorithm  delivers  rates  at 
Rfhx  =  1  Mbps.  The  average  performance  of  the  MODPC  algorithm  is  assessed  in  Figure  57. 


Rosslyn  at  1.9  GHz:  TX-2=4  m/s;  p  v=-149  dBW 
MODPC  Algorithm:  sc=10  dB;  /?=0.9;  h>mo=[0.99,0.01] 
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Figure  56.  MODPC  performance  for  mobile  TX-2  traveling  along  the  x  axis;  Rb,x  =  1  Mbps;  delay  r  =  0.9 
ms;  frame  duration  tF  =  1  ms;  fc  =  1-9  GHz;  coding  gain  gc  =  10  dB.  Upper  Panel:  Link  power  along  the 
track;  Lower  Panel:  Link  rates  along  the  track. 
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Figure  57  assesses  the  average  performance  of  the  MODPC  algorithm  as  the  mobile  network  travels 
along  the  x  axis.  Also  plotted  is  an  upper  bound  on  performance  given  by  the  Exact  FM  algorithm  (de¬ 
scription  forthcoming).  Neither  algorithm  meets  the  QoS  constraints  at  every  point  of  the  x-axis  track. 
Therefore,  the  next  two  paragraphs  detail  how  the  average  performances  account  for  this  loss  of  QoS. 
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Figure  57.  MODPC  performance  averaged  over  z-axis  track;  QoS  Rb,x  =  1  Mbps;  delay  t  =  0.9  ms;  frame 
duration  tF  =  1  ms;  fc  =  1.9  GHz;  coding  gain  gc  =  10  dB. 


MODPC  Algorithm:  let  Qmo  index  those  frames  meeting  the  QoS: 

Qmo  :=  {n  G  [1  :  NF]  -  1  :  7m(Pc;  tn)>lx]  tn  =  ntF}, 


where  NF  is  the  number  of  frames  along  the  track.  Averaging  the  rates  over  the  entire  track  mixes  the  un¬ 
successful  QoS  transmission  with  the  transmission  that  met  the  QoS.  This  mixing  is  problematic  when  one 
link  dominates.  Indeed,  Section  8.4  reports  Link  1  always  delivers  high  rates  although  Link  2  was  almost 
shut  down — total  throughput  is  high  while  the  QoS  is  never  satisfied.  Therefore,  the  average  network  rate 
is  defined  as  average  rate  only  when  the  network  delivers  QoS: 


Rmo 


1 

iVp 


E  iT 

»£Smo 


Rb^iin) 

Rb,2{tn) 


The  MODPC  algorithm  always  transmits  minimum  power  regardless  of  the  QoS.  Therefore,  network 
power  is  the  average  power: 

1  nf- i 

PMO  =  X]  lTPC'(^)- 

F  n= 0 


Exact  FM  Algorithm  :  This  algorithm  operates  with  zero  system  delay,  global  network  knowledge  (i.e.,  the 
gain  matrix),  and  solves  for  minimum  power  pmin  meeting  the  QoS  constraints.  Not  every  position  on  the 
.'/•-axis  track  admits  a  solution.  Let  Qxfm  index  those  frames  meeting  the  QoS: 

Qxfm  :=  {n  G  [1  :  Np]  -  1  :  7M(Pminfa))>7x;  tn  =  ntF}. 
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Define  the  average  network  rate  as  the  average  rate  only  when  the  network  delivers  the  QoS: 


-Rxfm  = 


1 

Nt 


£  11 


Rb,x 

Jh.x 


=  2  Rhx  x 


Qxfm 

The  average  power  used  by  the  exact  FM  algorithm  is  defined  as 

1 


Qxfm 

Nt 


PXFM  = 


I  Qxfm  | 


*52  lTPmin  (*n)- 


neQxFM 

Indeed,  the  power  is  not  defined  when  the  QoS  constraint  cannot  be  satisfied. 


Discussion:  Figure  57  reports  that  the  MODPC  algorithm  delivers  approximately  70%  of  the  best  possible 
average  network  rate  while  using  less  power.  The  MODPEC  algorithm  delivers  this  performance  against 
the  system  delay  of  r  =  0.9  ms  and  no  global  knowledge  of  the  network — each  link  uses  only  information 
measured  at  its  receiver.  The  next  section  assesses  the  MODPC  algorithm  against  larger  system  delays. 


10.2  MODPC  PERFORMANCE  AGAINST  DELAY 

Figure  58  compares  the  average  power-rate  performance  for  the  Exact  FM  algorithm  to  the  MODPC 
algorithm.  As  in  Section  10.1,  the  exact  FM  algorithm  operates  with  zero  system  delay,  global  network 
knowledge,  and  solves  for  minimum  power  pmin  that  meets  the  QoS  constraints.  Therefore,  its  perfor¬ 
mance  has  not  shifted  from  Figure  57.  The  MODPC  algorithm  “hangs”  tight  until  the  delay  exceeds  8 
ms.  At  9  ms  delay,  the  power  jumps  by  approximately  12  dB.  Increasing  delays  past  9  ms  further  increase 
the  power  while  throughput  rolls  off.  Thus,  the  MODPC  algorithm  has  a  delay  threshold — exceeding  this 
threshold  incurs  system  collapse. 
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Figure  58.  Average  MODPC  performance  along  z-axis  track  as  a  function  of  delay;  Rb,x  =  1  Mbps;  frame 
duration  tF  =  1  ms;  fc  =  1.9  GHz;  coding  gain  gc  =  10  dB. 
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11.  INTERLUDE:  KALMAN  FILTER 


A  typical  Kalman  filter  for  wireless  power  control  usually  tracks  the  SINR  and  update  power  based  on 
the  predicted  SINR.  This  section  reviews  the  Kalman  filter  estimating  the  SINR.  The  resutling  state  space 
is  used  in  the  forthcoming  sections. 

11.1  KALMAN  FILTER  MODELS 

The  simplest  discrete  Kalman  filter  is  based  on  the  state-space  model  and  measurement  model  of  the 
form  [5,  Equations  4.2.1-4.2.2]:  for  t  =  0,  1,  . . . 

Xt+1  =  <fhx4  +  wt,  (18) 

Z  t  =  Htxt  +  \t.  (19) 

The  transition  matrix  4?/  and  the  measurement  matrix  Ht  are  known.  The  process  noise  { wt }  and  the  mea¬ 
surement  noise  { V/ }  are  zero-mean,  white  noise  sequences  that  arc  uncorrelated  and  also  uncorrelated  with 
the  initial  state  [5,  Equations  4.2.3-4.2.5J 

Txol  1  HPo  o  o' 

E  w  t  [xj  w  J,  \J,]  =  0  Stjt'Qt  0  .  (20) 

_  Lv*  J  J  L°  °  St,t'Rt_ 

The  process  covariance  {Qi}  and  the  noise  covariance  { Rt }  arc  known.  The  Kalman  filter  advances  the 
prior  estimate  of  the  state  x~t  and  the  prior  estimate  of  its  error  covariance  Pf  from  the  current  measure¬ 


ment  [5,  Equations  4.2.17,  4.2.8,  4.2.22]: 

I<t  =  Pf7 Hj  (Htp-Hj  +  Rty'  (21) 

%  =  %  +  Kt  (zt  -  Htx7 )  (22) 

Pt  =  (. I-KtHt)Pt ~  (23) 

and  then  updates  for  the  next  recursion  [5,  Equations  4.2.23,  4.2.25]: 

Pt+i  =  *tPt$J  +  Qt  (24) 

x"+1  =  $txt.  (25) 


11.2  CHANNEL  GAIN  AND  INTERFERENCE 


The  power  control  algorithms  in  this  report  assume  perfect  SINR  measurements.  This  assumption  al¬ 
lows  algorithm  comparison  without  the  confounding  effects  of  measurement  errors  and  disrupting  trans¬ 
missions  (e.g.,  shutting  off  transmissions  to  measure  the  interference  is  impractical.)  In  a  wireless  power 
control  system,  the  receiver  sets  its  transmitter’s  power.  This  knowledge  of  the  transmitter’s  power  makes 
SINR  estimates  possible.  Received  power  on  link  l  =  1, . . . ,  L 

Pr,i( t)  =  9u(t)pr(t )  +  Pi,i(t)  +  PN,i(f )  [watts] 

is  a  standard  measurement  in  most  so  ft  ware-defined  radios  (SDR).  Define  the  state  vector 
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Assume  the  state-space  model  of  the  form 


xi(tn+i)  ■=  xt(tn )  +  A xi(tn). 

The  measurement  model  for  the  received  power  is 

PR,l(tn)  =  [pr,i(tn)  1]  x(f„)  +pN,i{tn). 

The  transition  matrix  and  the  measurement  matrix  are  known  because  the  transmitter’s  power  is  controlled 
by  the  receiver: 


m)  =  i 

Hl{t )  =  [pr,l(t)  l]  . 

Assume  the  additive  RF  noise  is  constant  for  all  receivers: 

Rl(tn)  =  cr2N. 


The  covariance  of  the  process  noise  is 


Qlitn)  =  E[Axi{tn)J  Axi{tn)} 


E[Agfi{tn)\ 

E[ApIti(tn)Agu(tn)] 


E[Agii{tn)Apn(tn)\ 

E[Ap2ij(tn )] 


Assume  there  is  no  delay  or  error  between  the  control  power  pcj  (t)  at  the  receiver  and  the  transmitted 
power: 

PT,i{t)  =  Pc,l(t )• 

As  the  Zth  link  increments  to  time  tn,  assume  prior  estimates  x;~(t)  and  Pt~  (t).  Update  these  prior  esti¬ 
mates  as  follows: 

Ki{tn )  =  Pi(tn)~ +  Ri(tn)j 
%{tn)  =  %{tn)~  +  Ki(tn)  ( PR,l(tn )  -  Hi(tn)%(tny) 

mn)  =  {I  -  K^tyHidnWitn)- . 


Update  or  predict  for  the  next  recursion  as  follows: 


Pl(tn+ 1)  =  Pi(tn)  +  Ql{tn) 

Xl(tn+ 1)~  = 


The  SINR  at  the  start  of  the  nth  frame  is  estimated  as 

/,  s  9ll(tn)PT,l(tn )  gii(tn)pT,i(tn) 

iMptn)  - TA\  I - ~  -  u  \  , - 2T- 

PI, l  ( tn  )  +  PN,l  {tn  )  Pi, l  ( tn  )  +  (?N 

Knowing  the  covariances  is  a  basic  problem  for  the  Kalman  filter.  The  next  sections  use  the  H2/H°°  ap¬ 
proach  to  drop  the  covariance  assumptions. 
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12.  H 2  ALGORITHMS 


Chaves,  Abbas-Turki,  Abou-Kandil,  and  Romano  approach  power  control  using  [13]:  “a  standard 
target  tracking  algorithm  with  a  dynamic  target  QoS”.  This  section  reviews  their  Opportunistic  QoS  ap¬ 
proaches  and  assesses  the  performance  of  their  static  H2  algorithm  (7A-DPC)  and  their  adaptive  H2  al¬ 
gorithm  (A-TQ-DPC)  on  the  two-link  mobile  wireless  scenario  of  Section  8.  The  FM  algorithm  updates 
power  on  the  l-  th  link  using  the  difference  between  the  external  SINR  7 x,i  that  is  fixed  throughout  the 
wireless  operation  and  the  measured  SINR  7 x,i{tn)  that  varies  in  the  mobile  scenario  [22]: 

Pc,l(tn+ 1)  =  Pc,l(tn)  +  (1  -  A)  {lx, l  -  lM,l(tn)}  ■ 

The  step-size  for  the  l- th  link  is  denoted  by  0  <  A  <  1-  The  dB  scaling  is  denoted  by  the  overline: 

fiefon)  =  10  log  10(pc,l(tn))- 

The  H 2  algorithms  update  power  as 

Pc,l(tn+ 1)  =  Pc, l  (A)  +  (1  -  A)  {lT,l{tn)  —  7m,z(A)}  ,  (26) 

where  the  target  SINR  7 Ti(fn)  varies  to  deliver  [13]:  “high  flexibility  on  the  quality-of-service  (QoS)  pro¬ 
vision  by  exploiting  link  quality  variation.”  The  variation  of  the  target  SINR  “steers”  state-space  model 
tracking  the  power  and  rate  along  a  trajectory  realizing  an  optimal  trade-off  between  power  and  rate.  More 
formally,  the  II2  algorithms  minimize  a  quadratic  form  of  the  power  and  QoS. 

This  section  starts  by  reviewing  the  state-space  model  and  the  quadratic  form.  Urban  simulations  arc 
run  using  nominal  values  that  parametrize  the  quadratic  form.  These  simulations  provide  a  concrete  ex¬ 
ample  to  discuss  the  quadratic  form.  This  section  concludes  by  discussing  parameterizations  the  quadratic 
form  to  deliver  better  performance.  A  state-space  model  determines  the  evolution  of  the  target  SINR  using 
a  state  vector  measuring  the  deviation  from  the  minimum  power  [13,  Equation  11] 

ef(A)  ■■  =  Pc,l(tn)  ~  Pm  in 

and  the  deviation  from  the  maximum  SINR  [13,  Equation  10] 


lT,l(tn )  7max  +  ul(tn), 

where  ui(tn )  is  the  control  signal.  Subtracting  pmin  from  Equation  (26)  determines  state  equation  of  the 
power  derivation 

ef(Wl)  =  A ef(tn)  +  (1  -  A)  [lT,l{tn)  -  filitn)  -  £>min}  , 

where  the  channel  gain  [13,  Equation  2] 


in(tn)  '■= 


_ 9ll(tn) _ 

PN,l(tn)  T  9lm{tn)PC,m{tn 


t) 


is  computed  from  the  measured  SINR  and  the  known  power.  Therefore,  the  state  vector  for  the  l- th  link 
[13,  Equation  12] 


X/0n) 


lT,l(tn)_ 


has  the  state  equation  [13,  Equation  13] 


Xz(tn+l)  =  Axfitn)  +  +  B2Ul(tn), 


73 


structured  by  the  matrices  [13,  Equation  14]: 


A  = 


'P  (1-/3) 

0  0 


0  -(1-/3)  -(1-/3)'. 

10  0 


B-2 


0 

1 


and  driven  by  the  exogenous  input  [13,  Equation  14]: 


w  i(tn) 


7max 

blifn) 

_  P min  . 


The  quadratic  trade-off  between  power  and  QoS  seeks  a  control  sequence  ui(tn )  to  minimize  the  asso¬ 
ciated  quadratic  form  [13,  Equation  15] 


N 

Q(u)  :=  ^xj (tn)Qxi(tn)  +  r(tn)ui(tn)2\  Q 

n= 0 


1  o 
0  0 


(27) 


determined  by  the  sequence  of  weighting  parameters  r(tn)  >  0.  Small  values  of  the  weighting  parame¬ 
ters  force  the  power  to  track  the  minimum  power  pnnn  while  accepting  a  lower  QoS.  Larger  values  of  the 
weighting  parameters  force  a  better  QoS  at  the  cost  of  more  power.  The  challenge  is  to  adapt  the  trade-off 
to  the  time- varying  gain  matrix  or,  equivalently,  an  RF  scenario.  Chaves  solves  the  quadratic  minimization 
by  converting  to  a  feedback  control  problem  [13,  Equation  16] 


q  .  f  *l(tn+ 1)  =  Axi(tn)  +  BiWl{tn)  +  B2Ul(tn ) 

T'\  Z  i(tn)  =  CiXi(tn)  +  D12ui(tn) 

that  is  structured  by  the  matrices  [13,  Equation  14]: 


Ci 


1  0 
0  0  ’ 


D 12 


0 

r(tn)1/2 


Figure  59  illustrates  the  feedback  control  determined  by  the  system  G.  Chaves  produces  the  static  con¬ 
troller  by  setting  r(tn)  =  r  constant  and  N  — >  oo  to  obtain  the  control  signal  [13,  Equation  18,  19] 

ui(tn)  =  - Kxi(tn )  =  -  jr  +  BJPB2 1  BjPA 

where  P  is  the  solution  to  the  discrete-time  algebraic  Riccati  equation  [13,  Equation  20] 

P  =  AJPA  +  Q  -  ATPB2{r  + B]PB2}  1  BjPA. 


Jf(fn) 
7  Ti(^n) 


Algorithm  3  {U2- DPC)  Select  the  step  size  0  <  f5  <  1  and  QoS  weight  r  >  0.  Compute  the  solution  P  to 
the  discrete -time  algebraic  Riccati  equation  [13,  Equation  20] 

P  =  AJPA  +  Q  -  ATPB2{r  + B]PB2}  1  BjPA. 

Compute  the  feedback  gain  matrix  [13,  Equation  19] 

I\  =  \r  +  Bj  PB2  }  1  Bl  PA. 
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>z 


Figure  59.  Kalman  feedback  for  QoS  power  control. 


'  4(t o) ' 

Pini,/  Pmin 

T  max 

Select  the  minimum  power  pmin,  the  initial  transmission  power  pini  j,  and  the  maximum  SINR  7max.  Initial¬ 
ize  the  state  vectors 

xl(to)  = 

on  all  the  links  l  =  1,  . . . ,  L.  Assume  the  wireless  network  has  been  operating  with  initial  power  up  to  the 
time  the  power  control  is  turned  on:  pc,i(t )  =  Pini, l  for  t  <  to  and  for  l  =  1,  . . . ,  L 

H-l  The  network  is  synchronized  and  all  updates  occur  at  the  start  of  each  frame:  tn  :=  ntF- 
H-2  Each  receiver  measures  its  SINR  at  the  start  of  each  frame: 

9ll(fn)PC,l(fn  -  t) 


' YM,l(tn )  = 


PN,l(fn)  ~F  y^ivvyi  Ql.m.(tn)PC.m.(tn  '7~) 


H-3a  Extract  the  adaptive  target  SINR  from  the  state  vector 

'  *?(<«: ' 


X/(*n)  = 


lT,l{tn)- 


H-3b  Update  the  exogenous  input  vector  from  the  measured  SINR 


1M,1  (tn 


7max 

Tl(tn) 
.  .Pmin 


=  w  i(tr, 


H-3c  Update  the  state  vector 

xi(tn+i)  =  (A  +  B2K)xi(tn)  +  Bpwftn). 
H-4  The  power  is  updated  for  the  frame  [tn,  tn+ 1)  as 


PC,l(tn+ 1) 


PC,l(tn ) 


7 T,l(tn)  \  1  13 
7 M,l(tn)  J 
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12.1  U2- DPC  MOBILE  PERFORMANCE  AT  1.9  GHZ 


Table  1 1  lists  the  RF  parameters  of  the  gain  matrix,  the  SINR,  and  the  updating  of  the  wireless  control 
algorithms.  These  parameters  are  outside  the  control  of  the  T~L2- DPC  algorithm.  Table  14  lists  the  inputs 
to  the  772-DPC  algorithm.  The  QoS  is  1  Mbps  on  both  links;  the  large  weighting  parameter  r  biases  the 
power-rate  trade-off  to  the  QoS;  the  frame  duration  of  tp  =  1  ms  sets  the  power  control  update  rate;  the 
delay  is  a  fraction  of  a  frame. 


Table  14.  %2-DPC  Parameters. 


Variable 

Value 

Description 

Rb,X 

Pmin 

P  max 

Pini 

P 

r 

1  Mbps 

0  W 

1  W 

1 0  //W 
0.9 

500 

Throughput 

Minimum  power 
Maximum  power 

Initial  power 

Stepsize 

Weighting  parameter 

Figure  60  reports  the  power  and  the  rate  for  the  network  during  the  time  TX-2  travels  along  the  x-axis 
track.  The  upper  panel  compares  the  power  on  each  link  computed  by  the  distributed  Ti2- DPC  algorithm. 
The  lower  panel  compares  throughput  on  each  link.  Although  the  coding  gain  is  gc  =  10  dB,  throughput 
is  set  using  the  capacity:  Rb(t)  =  f  n  x  log2(l  +  7m (!))■  On  the  average,  the  T/J-DPC  algorithm  delivers 
rates  at  !h,.x  =  1  Mbps.  The  next  plot  reports  the  average  performance  along  the  x-axis  track. 


Figure  60.  "H2-DPC  performance  along  x-axis  track;  Rt,x  =  1  Mbps;  frame  duration  tF  =  1  ms;  fc  =  1-9 
GHz;  coding  gain  gc  =  10  dB.  Upper  Panel:  power  on  each  link;  Lower  Panel:  throughput  on  each  link. 
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Figure  61  assesses  the  average  performance  of  the  PI2- DPC  algorithm  as  the  mobile  network  travels 
along  the  x-axis  track.  As  described  in  Section  9. 1,  an  upper  bound  on  performance  given  by  the  Exact  FM 
algorithm  is  plotted.  Neither  algorithm  can  meet  the  QoS  constraints  at  every  point  along  the  x-axis  track. 
The  next  paragraph  details  how  the  average  performances  of  both  algorithms  account  for  this  loss  of  QoS. 
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Figure  61.  Average  PL2- DPC  performance  along  x-axis  track  with  delay  r  =  0.9  ms;  RbtX  =  1  Mbps;  frame 
duration  tF  =  1  ms;  fc  =  1.9  GHz;  coding  gain  gc  =  10  dB. 


PL2 -DPC  Algorithm:  Let  Qy2  denote  those  frames  that  meet  the  QoS 

Qn 2  :=  {ne[l:  Nf\  -  1  :  7Af(Pc(*n);in)>7x;  tn  =  ntp}, 


where  Np  denotes  the  number  of  frames  along  the  track.  Averaging  the  rates  over  the  entire  track  mixes 
the  unsuccessful  QoS  transmission  with  the  transmission  that  met  the  QoS.  This  mixing  is  problematic 
when  one  link  dominates.  For  example,  Section  8.4  reports  Link  1  delivers  high  rates  while  Link  2  is  al¬ 
most  shut  down.  Consequently,  total  throughput  is  high  although  QoS  is  never  satisfied.  Therefore,  the 
average  network  rate  will  be  defined  as  average  rate  only  when  the  network  delivers  QoS: 


H-H'2 


1 


E  iT 


n£QH2 


R-b,l  (fn) 
Rb,  2  (t?i) 


Average  network  power  will  be  defined  as  the  average  total  power 


1  Np  —  l 

Pn2  =  J2  lTP c(tn) 

F  n= 0 

because  the  Tf2-DPC  algorithm  always  transmits  a  minimum  power. 
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Exact  FM  Algorithm  (repeated  for  convenience):  The  exact  FM  algorithm  operates  with  zero  system  de¬ 
lay,  global  network  knowledge  of  the  gain  matrix,  and  power  levels  to  solve  for  minimum  power  pmin  that 
meets  the  QoS  constraints.  Not  every  position  along  the  rc-axis  track  admits  a  solution.  For  the  exact  FM 
algorithm,  let  Qxfm  index  those  frames  meeting  the  QoS 

Qxfm  :=  {n  G  [1  :  Np]  -  1  :  7M(pmmfa);  tn)>Ex]  tn  =  nip}. 


Define  the  average  network  rate  as  the  average  rate  only  when  the  network  delivers  the  QoS: 


Rxfm  =  xr 

JV t 


E 

«SQxfm 


Rb,X 

Rb.x 


=  2  Rhx  x 


Qxfm 
Nt  ' 


Therefore,  the  average  power  used  by  the  exact  FM  algorithm  is 


PXFM  = 


1 

I  Qxfm  | 


lTPmin(*n). 

neQxFM 


Discussion:  Figure  61  reports  that  the  "H2-DPC  algorithm  delivers  approximately  half  of  the  best  possible 
network  rate  using  only  slightly  more  power.  The  "H2-DPC  algorithm  delivers  this  performance  against  the 
system  delay  of  r  =  0.9  ms  and  no  global  knowledge  of  the  network — each  link  uses  only  information 
measured  at  its  receiver.  The  next  section  assesses  the  'H2-DPC  algorithm  against  larger  system  delays. 

12.2  n2- DPC  PERFORMANCE  AGAINST  DELAY 

Figure  62  compares  the  average  power-rate  performance  for  the  Exact  FM  algorithm  to  the  "H2-DPC 
algorithm.  As  in  Section  12.1,  the  exact  FM  algorithm  operates  with  zero  system  delay,  global  network 
knowledge,  and  solves  for  minimum  power  pmin  meeting  the  QoS  constraints.  Therefore,  its  performance 
has  not  shifted  from  Figure  61.  The  "H2-DPC  algorithm  “hangs”  tight  until  the  delay  exceeds  8  ms.  At  9 
ms  delay,  the  power  jumps  by  25  dB.  Increasing  delays  in  excess  of  9  ms  force  the  power  to  increase  while 
the  throughput  rolls  off.  Thus,  the  EC2 -DPC  algorithm  exhibits  delay  threshold — exceeding  this  threshold 
incurs  system  collapse. 

12.3  ADAPTIVE  EL2-DPC 

Chaves  points  out  that  adapting  the  weighting  parameters  r{tn)  >  0  to  the  time-varying  channel  pro¬ 
duces  an  adaptive  "H2-DPC  algorithm  and  offers  the  following  rule  as  a  starting  point  [13,  Equation  21]: 

r(tn)  ■■=  max-0  x  (7max  -  rmin},  (28) 

where  rm;n  >  0  and  <j)  >  0  selects  throughput  or  power  savings  via  <i>  “large”  or  “small,”  respectively.  Min¬ 
imizing  the  quadratic  for  this  adaptive  r(tn)  forces  a  time- varying  gain  K(tn )  computed  from  a  backwards 
Riccati  [13,  Equation  22] 

I\(tn)  =  (r(tn)  +  B2P{tn+1)B2)~1BjP(tn+i)A 

and  a  backwards  Riccati  for  the  projection  matrix  P{tn).  These  backwards  recursions  are  not  possible  be¬ 
cause  the  weighting  parameters  r(tn)  arc  computed  from  propagation  conditions.  Chaves  works  around 
this  future -past  conflict  as  follows  [13]:  “Therefore,  the  adaptive  EL2- DPC  algorithm  consists  in  imple¬ 
menting  the  EL2- DPC  algorithm  for  a  different  r(tn)  for  each  time  slot,  according  to  Equation  (28).  We  call 
this  procedure  the  adaptive  EL2  distributed  power  control  (A-"H2-DPC)  algorithm.” 
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Rosslyn  at  1.9  GHz:  TX-2=4  m/s;  pN=- 149  dB 


Figure  62.  Average  H2- DPC  power-rate  performance  along  i-axis  track  subjected  to  delays;  Rb.x  =  1 
Mbps;  frame  duration  tF  =  1  ms;  fc  =  1.9  GHz;  coding  gain  gc  =  10  dB. 


Algorithm  4  (A-772-DPC)  Select  the  step  size  0  <  f3  <  1.  Select  rmjn  >  0.  Select  d>  >  0  “large  ” 
for  throughput  or  <fi  >  0  “small”  to  save  power.  Select  the  minimum  power  pmin,  the  initial  transmission 
power  p]u]j,  and  the  maximum  SINR  7max-  Initialize  the  state  vectors 

Pini,l  —  iPmin 

7  max 

on  all  the  links  l  =  1 ,  . . . ,  L.  Assume  the  wireless  network  has  been  operating  with  initial  power  up  to  the 
time  the  power  control  is  turned  on:  pc,i(t )  =  Pini,lfor  t  <  to  and  for  l  =  1,  . . . ,  L: 


xi(to)  = 


5(*>) 

It, o) 


H-l  The  network  is  synchronized  and  all  updates  occur  at  the  start  of  each  frame:  tn  :=  nt  p . 

H-2  Each  receiver  measures  its  SINR  at  the  start  of  each  frame: 

,  S  =  _ 9ll{tn)pc,l{tn  -  r) _ 

7M ,l  n  PN,l(tn )  +  gim(tn)pc,m(tn  ~  r)  ' 

H-3a  Extract  the  adaptive  target  SINR  from  the  state  vector 


xl(tn )  = 


7  T,l(tr 


lT,l{tn )• 


H-3b  Update  the  exogenous  input  vector  from  the  measured  SINR 


7max 

Pl{tn) 

.Pmin 


=  W  l(tn)- 
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H-3c  Compute  the  weighting  parameter  [13,  Equation  21  ]: 


ri(tn)  =  max{0  X  (7 max,/  -  rmin}- 

H-3d  Compute  the  solution  Pi(tn )  to  [13,  Equation  20] 

Pi(tn)  =  Pi(tn)A  +  Q  —  ATPi(tn)B2  |n(in)  +  BjPi(tn)B2 j*  Bj Pi(tn)A. 

H-3e  Compute  the  feedback  gain  matrix  [13,  Equation  19] 

Klitn )  =  |  ri(tn)  +  BjPi(tn)B2 1  bJ  Pi{tn)A. 


H-3f  Update  the  state  vector 


xi{tn+1)  =  (A  +  B2Kl(tn))xi(tn )  +  Biw(tn). 

H-4  The  power  is  updated  for  the  frame  [tn ,  tn+\ ')  as 

u  ^  u  ^  f  lT,l{tn)  )1_/3 

PC,l{tn+ 1)  =  VC,l(tn )  <  - 7~T  ) 

Figure  63  shows  the  adaptive  TL2  performance  is  comparable  to  the  non-adaptive  algorithm  under  moder¬ 
ate  delay.  Figure  64  reports  that  the  the  adaptive  TL2  breaks  down  with  the  delay  exceeds  8  ms.  Figure  65 
is  the  power  allocation  of  this  adaptive  algorithm.  This  figure  shows  that  the  excessive  delay  causes  this 
algorithm  to  “chatter”  and  instability  is  likely  held  in  check  by  the  power  limiter. 


Rosslj  n  at  1.9  GHz:  TX-2=4  m/s;  p  =-149  dB 
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Figure  63.  Average  A-7f2-DPC  performance  along  x-axis  track:  r  =  0.9  ms. 
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Rosslvn  at  1.9  GHz:  TX-2=4  m/s;  p  =- 149  dB 


Average  Network  Power  (dBW):  Rb  1  Mbps;  t=9  ms;  1^=1000  /is 


Figure  64.  Average  A-%2-DPC  performance  along  i-axis  track:  r  =  9  ms. 
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Figure  65.  Power  allocation  for  the  A-%2-DPC  along  x-axis  track:  r  =  9  ms. 
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13.  THE  FM-TD  ALGORITHM 


The  preceding  sections  show  the  limitations  of  the  FM  algorithm,  Kalman  filtering,  the  H 2  and  H°° 
approaches  that  model  the  errors  caused  by  the  time  delay  as  noise.  Increasing  time  delays  causes  all  these 
approaches  to  fail — time  delay  is  not  part  of  the  modeling.  Gunnarsson,  Guustafsson,  and  Blom  produced 
a  superb  analysis  of  time-delay  effects  [29] : 


“. . .  round-trip  delays  in  the  power  control  loop  result  in  oscillations.  Essentially,  the  core 
problem  is  that  the  measurements  do  not  reflect  the  most  recent  power  updates.  However, 
these  arc  known  to  the  algorithm,  and  can  be  compensated  for.  Therefore,  we  adjust  the  mea¬ 
sured  SIR  according  to  the  power  control  commands  that  have  been  sent  but  whose  effects 
have  not  yet  been  experienced  by  the  receiver.” 


The  problem  is  that  each  receiver  measures  the  time-delayed  SINR  at  the  start  of  each  frame: 

(  ,  =  _ gnjtnjpcAtn  ~  T ) _ 

7A/,i  n  PN,l(tn )  +  Em/1  9lm(tn)PC,m(tn  ~  r) 

whereas  the  update  should  use  the  SINR  without  time  delay 

9ll(tn)Pc,l(tn ) 


lM,i{tn)  = 


PN,l(tn)  9lm(tn)PC.m(tn 


The  solution  is  approximate  7 M,i(tn)  using  the  measured  SINR  and  the  past  power  control  pc(t„  —  r) 
stored  at  the  receiver  [29] : 


7 


tn ) 


9ll(tn)pc,l{tr 


PN,l{tn)  +  Zm*l9lm(tn)pC,m(tn  T) 
PC,l(tn ) 


=  7 M,l(tn) 


PC, l  (tn  -  t)‘ 


Thus,  the  FM  algorithm — and  the  other  wireless  power  control  algorithms — only  needs  a  slight  modifica¬ 
tion  to  compensate  for  time  delays. 


82 


Algorithm  5  (FM-TD)  [29]  All  receivers  have  an  estimate  Tp  of  the  system  delay  t. 


FMTD-1  The  network  is  synchronized  and  all  updates  occur  at  the  start  of  each  frame:  tn  :=  nt  p . 

FMTD-2  Each  receiver  measures  the  time-delay  SINR  at  the  start  of  each  frame: 

,  S  =  _ 9ll{tn)PC,l{tn  -  t) _ 

lM'1  n  PN,l(fn )  +  Em±Z  9lm(tn)pc,m(tn  ~  t)  ' 

FMTD-3  Each  receiver  compensates  for  the  time  delay  in  the  SINR: 

PC,l(fn) 


lM,l{jF\  tn)  =  2/M,l(tr 


PC, l  (fn  ~TF)' 


FMTD-4  The  power  is  updated  for  the  frame  [tn,  tn+\)  using  the  delay-compenstated  SINR: 

1-/3 


PC,l{tn+ 1)  =  PC, l  (fn) 


lX,l 


lM,l(fF',tn ) 


FMTD-5  The  control  power  is  limited 

PC,lifn+l)  =  [Pmin  ±  PC,l(tn+ 1)  ^  Pmax] • 


Figure  66  reports  this  delay  compensation  for  the  two-link  mobile  wireless  scenario.  Instead  of  the 
power-rate  performance  rolling  off  with  increasing  delay  as  seen  in  the  FM  algorithm  (Figure  55)  and  the 
H2  algorithm  (Figure  62),  the  FM-TD  does  compensate  for  time  delays  sweeping  out  to  20  ms.  However, 
this  performance  assumes  the  system  delay  is  known  exactly: 


T  =  Tp. 


The  natural  question  asks:  How  does  FM-TD  performance  degrade  when  the  system  delay  r  is  not  known 
exactly? 

Both  the  FM  algorithm  and  the  H2  algorithm  breakdown  when  the  system  delay  exceeds  r  =  10  ms 
(Figure  55  and  Figure  62).  Accordingly,  the  system  delay  fixed  at  r  =  10  ms  and  Tp  is  vailed  around  this 
delay: 

Tp  =  T  x  (1  ±  25%)  G  [7.5, 12.5]  [ms]. 

Figure  67  reports  that  the  FM-TD  is  robust  with  respect  to  errors  in  the  delay  of  ±25%:  Although  30  Tp’ s 
were  generated,  rounding  to  the  frame  increments  produced  the  discrete  distribution  of  FM-TD  power-rate 
performance. 
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Rosslvn  at  1.9  GHz:  TX-2=4  m/s;  p  =- 149  dB 


Average  Network  Power  (dBW):  Rh  >  1  Mbps;  tf=1000  g$ 


Figure  66.  Average  FM-TD  performance  along  x-axis  track  subjected  to  delays;  QoS  RbjX  =  1  Mbps; 
frame  duration  tF  =  1  ms;  fc  =  1-9  GHz;  coding  gain  gc  =  10  dB;  delay  is  known:  r  =  tf. 


1.9 


1.8 


??1.7 


5  1.6 
ck 

* 

o 

I  1.5 


1  1.4  - 


C 


1.3 


1.2 
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Figure  67.  Average  FM-TD  performance  along  x-axis  track  subjected  to  delays;  QoS  Rh.x  =  1  Mbps; 
frame  duration  tF  =  1  ms;  fc  =  1-9  GHz;  coding  gain  gc  =  10  dB;  delay  is  approximate:  tf  «  r. 
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14.  DELAY  AND  STABILITY 


The  preceding  sections  show  increasing  system  delay  causes  the  algorithms  to  break  down.  This  break¬ 
down  is  analyzed  by  considering 

•  Stability  of  the  FM  algorithm 

•  Variation  of  the  channel 

•  Effect  of  delay  on  stability. 

The  FM  algorithm,  without  delay  and  step  size  is  [18,  Equation  6]: 

p(4+l)  =p(*n)  y 

7 M\tn ) 

Assuming  the  channel  does  not  change,  this  distributed  recursion  is  equivalent  to  the  power  method  [18, 
Equation  5] 

p(fn+i)  =y  +  Fp(fn),  (29) 

where  F  is  the  feedback  matrix  F  =  G^'YxGj  determined  by  the  channel  gain  matrices  and  Y\-  := 
diag(7x)  Assuming  spectral  radius  of  the  feedback  matrix  is  strictly  less  than  unity  p(F)  <  1,  the  FM 
algorithm  converges  to  the  solution  of  (/  —  F)p  =  y  by  the  partial  sums 

n 

p(fn)  =  ^Ffcy^(/-F)-1y. 

k= 0 

Thus,  convergence  of  the  FM  algorithm  is  guaranteed  assuming  (/)  that  the  channel  does  not  change  and 
(ii)  the  spectral  radius  p{F)  <  1.  Figure  68  plots  the  spectral  radius  p(F)  along  the  re -ax  is  track.  This  plot 
shows  p(F)  <  1  excepting  a  few  points  in  the  initial  part  of  the  track  (0  <  x  <  100).  Thus,  the  FM  al¬ 
gorithm  converges  at  almost  all  points  along  the  track — provided  the  channel  does  not  change.  Therefore, 
the  discussion  turns  to  the  question:  What  effect  does  the  channel  variation  have  on  the  FM  convergence? 
The  absolute  error  between  the  matrix  inverse  (I  —  F)~ 1  and  the  associated  partial  sums  is  bounded  by 
Fenima  4: 

«FM  II  /r||nFM  +  l 

AbsErr(F;  nFM)  :=  (/  —  F)~l  —  ^  Fk  <  "  "  ;  n  =  0,l,.... 

k= 0 

Define  the  relative  error  between  the  matrix  inverse  (I  —  F)_1  and  the  associated  part ial  sums  as  follows: 

1T?  ^  AbsErr(F;  nFM) 

RelErr(F; nFM)  :=  • 

Figure  69  estimates  the  number  of  iterations  nppj  for  the  FM  algorithm  to  converge  to  “2-digit”  accuracy: 

RelErr(F;  nFM)  <  10-2, 

using  the  standard  matrix  norm.  This  plot  report  that  convergence  requires  less  than  200  iterations  over 
most  of  the  slow  fading  (0  <  x  <  100)  paid  of  the  track.  Over  the  fast  fading  part  of  the  track  (100  <  x), 
the  plot  reports  five  iterations  or  less  arc  needed. 
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Rosslyn  at  1.9  GHz:  Feedback  Matrix 


R  =  1  Mbps;  g  =10  dB;  p=- 149  dB 
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Figure  68.  Spectral  radius  of  the  feedback  matrix  a--axis  track;  QoS  Rb>x  =  1  Mbps;  fc  =  1.9  GHz;  coding 
gain  gc  =  10  dB. 
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Figure  69.  FM  algorithm  convergence  along  ;r-axis  track;  QoS  R.b,x  =  1  Mbps;  fc  =  1.9  GHz;  coding  gain 
gc  =  10  dB;  “blanks”  mark  where  p{F(x))  >  1. 


86 


Section  8  defines  the  channel  correlation  function  |C'^/l(Ax)|  and  the  correlation  distance  Ax (90%)  as 
the  biggest  distance  where 

\Chh(&x)\  >  0.9 

for  all  Ax  6  [0,  Ax(90%)].  Table  15  reports  this  correlation  distance  in  Column  2  for  both  the  slow-fading 
region  (0  <  x  <  100)  and  the  fast-fading  region  (x  >  100).  Column  3  lists  the  time  for  transmitter  TX-2 
to  decorrelate  when  traveling  at  v  =  4  m/s.  The  FM  algorithm  is  updating  at  the  frame  rate  of  t p  =  1 
ms.  Columns  4  and  5  lists  the  time  and  number  of  iterations  for  for  the  FM  algorithm  to  deliver  “two-digit” 
accuracy.  Table  15  and  shows  the  FM  algorithm  will  converge  to  “two-digit”  accuracy  in  both  the  slow- 
fading  and  fast-fading  regions — without  time  delays. 

Table  15.  Spatial  correlation  and  convergence  along  the  x-axis  track;  Transmitter  TX-2  traveling  at  v  =  4 
m/s;  frame  duration  tF  =  1  ms. 
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Figure  70  shows  the  FM  algorithm  performance  along  the  track  with  the  system  delay  of  r  =  10  ms. 
This  figure  shows  that  the  FM  algorithm  breaks  down  in  the  slow-fading  region  but  recovers  to  deliver  the 
QoS  in  the  fast-fading  region. 


Rosslyn  at  1.9  GHz:  TX-2=4  m/s;  />  v=-149  dBW 
FM  Algorithm:  Coding  Gain=10  dB;  /3= 0.9 
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Figure  70.  FM  algorithm  against  delay  r  =  10  ms  along  the  x-axis  track;  QoS  Rb,x  =  1  Mbps;  fc  =  1-9 
GHz;  coding  gain  gc  =  10  dB. 
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Time  delay  induces  instabilities.  The  FM  algorithm  without  stepsize  operating  in  a  static  channel 


Pl(tn+ 1)  =  Pl(tn ) 


lX,l 


lM,l(tn ) 
Pl(tn)  IX, l 


Pi(tn  -  r)  gu 


PN ,1  +  ^  9mlPm(tr 

rriy^l 


-  T 


(30) 


When  the  time  delay  is  zero  r  =  0,  the  leading  term  is  unity  and  this  formulation  reduces  to  Equation  (29): 


Pl(tn+ 1)  =  — ~PN,l  +  9mlPm(tn ) 

/T/7  '  fill 


IX, l 


gu 


m^l 


gu 


Analysis  of  the  single-link  system  suffices  to  show  how  delay  affects  the  system  stability.  Figure  7 1  is  a 
block  diagram  of  the  single-link  system. 


The  FM  algorithm  for  a  single  link  with  delay  and  stepsize  is 

T)/fJ 

where  g  is  the  the  channel  magnitude.  Converting  to  log  scale  and  approximating  r  ~  tj?k  produces  the 
linear  recursion 

p(n  +  1)  =  p(n )  -  (1  -  P)p{n  -  k)  +  y. 

This  recursion  is  stable  provided  the  roots  of 

ip(k,  P;  z )  =  zk+1  -  zk  +  (1  -  $)  =  0 

are  contained  in  the  closed  unit  disk  \z\  <  1  and  any  root  on  the  boundary  of  the  disk  \z\  =  1  is  simple  [71, 
Theorem  7.2.9. 6].  For  each  delay  k,  let  zm(k ,  /3)  denote  the  roots  of  the  recursion  m  =  1, . . . ,  (k  +  1): 

ijj(k,P;zm(k,P))  =  0. 

Figure  72  plots  the  maximum  of  the  magnitude  roots  as  a  function  of  the  stepsize. 

P  ma,yi{\zm(k,  P)\  :  m  =  1, . . . ,  (k  +  1)}. 

This  plot  shows  that  increasing  delays  require  P  to  migrate  toward  1.  The  two-link  simulations  in  this  re¬ 
port  all  use  P  =  0.9.  Figure  72  reports  the  FM  algorithm  becomes  unstable  for  k  =  16  at  P  =  0.9.  All 


simulations  in  this  report  set  the  update  increment  at  ip  =  1  ms.  Thus,  the  single  links  become  unstable  at 
the  delay  r  =  16  ms.  The  two-link  networks  in  this  report  all  reported  instability  when  the  delay  r  >  8 
ms.  These  two-link  networks  are  two  single-links  coupled  by  then-  mutual  interference.  If  each  single  link 
becomes  unstable  at  t  =  16  ms,  instability  for  t  <  16  should  be  expected  when  these  unstable  links  are 
coupled.  Therefore,  stabilization  of  each  single  link,  such  the  time-delay  compensation  algorithm  (Sec¬ 
tion  13),  is  a  necessary  condition  for  network  stability. 


0 

Figure  72.  Stability  of  the  recursion  for  the  single-link  system;  maximum  of  roots  of  zk+1  -  zk  +  (1  -  fj)  =  0 
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15.  EXTENDING  WIRELESS  POWER  CONTROL 

Performance  of  the  wireless  power  control  algorithms  was  determined  from  the  three-dimensional  ur¬ 
ban  simulation  operating  under  the  assumptions  of  ideal  measurements  and  perfect  timing.  As  such,  algo¬ 
rithm  performance  was  assessed  against  channel  variations  and  system  delays  without  the  entanglements 
of  modeling  non-ideal  measurements.  This  final  section  proposes  a  similar  approach  to  encompass  non¬ 
ideal  measurements  and  multivariate  wireless  power  control  algorithms: 

•  Non-ideal  measurements 

•  Joint  optimization  and  cross-layer  optimization 

•  Wideband  power  control. 

These  simulations  propose  to  follow  the  same  pattern  as  the  simulations  in  this  report:  upper  bounds  on 
algorithm  performance  arc  obtained  by  modeling  “best  possible”  assumptions;  relative  performance  shifts 
between  contending  algorithms  ranks  of  the  algorithms. 

15.1  NON-IDEAL  MEASUREMENTS 

Assessment  of  the  wireless  power  control  algorithms  assumed  ideal  measurements:  perfect  SINR  mea¬ 
surements,  exact  power  control  commands,  and  perfect  timing  for  the  updates.  There  are  two  payoffs  for 
this  approach:  (/)  an  upper  bound  on  algorithm  performance  is  established  assuming  ideal  measurements, 
(ii)  algorithm  sensitivity  to  non-ideal  measurements  is  computable.  The  latter  item  requires  error  model¬ 
ing  encompassing  measurement  systems  currently  operating  in  specific  radio  systems.  A  remarkable  paper 
by  Jian,  Sidiropoulos,  and  Giannakis  undertakes  such  modeling  [36].  The  paper  start  by  quantifying  the 
necessity  of  accurate  measurements  for  wireless  power  control  [36] : 

•  “A  1-dB  reduction  in  local  power  estimation  error  enables  accommodating  five  more  users  in  a  sys¬ 
tem  for  a  fixed  outage  probability,  under  a  variety  of  power  control  schemes.”  [74] 

•  Adapting  “signal  constellations  according  to  both  the  instantaneous  receiver  signal  power  and  local 
mean  (shadow)  power  . .  .require  accurate  shadow  power  estimation  and  prediction.”  [19] 

The  paper  identifies  limitations  of  commercial  power  estimation  techniques  [36]: 

•  “A  simple  type  of  window-based  estimator,  namely  weighted  sample  average  estimators  of  local 
mean  power,  is  currently  employed  in  many  commercial  communication  systems  such  as  GSM,  and 
the  Motorola  Personal  Access  Communication  System  (PACS)  ”  [33],  [59] 

•  “Both  analysis  . . .  and  experiment  shows  that  the  mean-square  error  (MSE)  performance  of  these 
window  based  estimators  deteriorates  severely  when  the  window  size  increases  beyond  certain 
value.” 

The  paper  produces  a  Kalman  filter  for  estimating  local  channel  gain  [36] : 

“that  always  meets  or  exceeds  the  performance  of  window-based  estimators  . . .  and  provides  a 
host  of  other  side  benefits  as  well.” 

This  approach  influenced  other  papers  [68],  [72].  Therefore,  these  papers  provide  a  natural  approach  to 
rank  wireless  power  control  algorithms  operating  from  non-ideal  measurements. 
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15.2  CROSS-LAYER  OPTIMIZATION 


The  wireless  power  control  algorithms  considered  in  this  report  were  limited  only  controlling  the 
power  in  the  PHY-layer.  A  massive  literature  exists  on  boosting  wireless  power  control  by  jointly  opti¬ 
mizing  other  PHY-layer  controllers  or  working  across  the  MAC  layer.  Two  surveys  on  this  Cross-Layer 
optimization  are  [50],  [58]  while  the  specific  algorithms  are  listed: 

•  Power  Control  &  Beamforming:  [30],  [61] 

•  Power  Control  &  Channel  Selection:  [54],  [68] 

•  Power  Control,  Scheduling  &  Routing:  [44],  [67],  [51],  [57]. 

The  payoff  from  these  extra  “knobs”  in  controller  is  expected  to  be  substantial.  For  example,  beam 
steering — essentially  managing  the  network’s  inference — nearly  doubles  network  throughput  in  this  urban 
simulation  [16].  However,  enlarging  the  number  of  “knobs”  in  the  domain  of  the  optimizer  typically  in¬ 
creases  system  delay  and  decreases  system  stability.  Algorithms  that  encompass  the  MAC  layer  are  depen¬ 
dent  on  the  network  messaging.  This  messaging  may  not  be  supported  in  the  distributed  or  non-centralized 
wireless  power  control  for  tactical  networks.  Therefore,  before  implementing  joint  optimization  scheme 
that  operate  across  multiple  layers,  simulations  can  assess  both  the  performance  and  the  stability  under  the 
assumption  of  perfect  cross-layer  communication.  If  the  proposed  algorithms  deliver  significant  perfor¬ 
mance  shifts,  consideration  of  the  messaging  complexity  is  justified. 

15.3  BANDWIDTH 

The  wireless  algorithms  considered  in  report  are  narrowband.  Power  control  algorithms  operating  over 
wide  bandwidths  typically  employ  OFDM  waveforms  [69],  [42],  [78].  Local  geometry  and  the  wavelength 
determine  the  variations  across  the  band.  For  example,  Figure  73  plots  the  wideband  gain  matrix  (left)  and 
the  SINR  (right)  over  the  20-MHz  band  centered  at  1.9  GHz  when  the  mobile  transmitter  TX-2  is  located 
at  x  =  0. 


Figure  73.  Gain  matrix  at  1 .9  GHz  (left);  SINR  at  1 .9  GHz  (right). 
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The  gains  arc  relatively  flat  in  comparison  to  the  interference  generated  by  Link  1.  That  interference 
has  15-dB  swings  forcing  corresponding  variations  in  the  SINR.  The  SINR  is  computed  assuming  both 
transmitters  each  broadcast  a  total  of  pr  =  1  watt  over  the  entire  frequency  band  fj>  =  20  MHz  producing 
a  signal  spectral  density  vector 


Pr(f)  = 


—  I"1" 

fB  k 


[watts  •  Hz  ’]. 


This  signal  spectral  density  is  reported  as  pr  =  —13  dB/MHz  on  the  SINR  plot.  The  total  noise  noise 
power  is  pn  =  —158  dBW  and  arises  from  the  noise  spectral  density  of  pn  =  —149  dB/MHz.  The  SINR 
is  the  ratio  of  the  local  signal  to  noise  and  interference: 


7m(pt;/)  =  ■ 


Gs(f)PT(f) 


PnU)  +  Gi(f)pr(f)  ’ 


where  the  signal  and  interference  matrices  arc 


Gs(f) 


9n(f)  0 

0  922  (/) 


and  Gi(f) 


0  9i2(/) 

92i(f)  0 


These  plots  reveals  several  features:  first,  the  SINR  of  the  links  arc  separated  by  approximately  15  dB. 
Therefore,  a  10-dB  coding  gain  delivered  by  CDMA  is  one  option  to  close  this  gap.  Second,  the  channels 
are  frequency  selective  so  the  OFDM  is  a  candidate  waveform.  Kuehn,  Dekorsy,  and  Kammeyer  observe 
[43]: 


“  OFDM  (Orthogonal  Frequency  Division  Multiplex)  is  used  to  combat  the  frequency  selec¬ 
tivity  of  the  mobile  radio  channel.  Therefore,  each  subcarrier  is  affected  by  flat  fading,  hence  a 
one-tap  equalizer  suffices  for  eliminating  channel  distortion.” 

Gueguen,  Crussiere,  and  Helard  observe  that  combining  OFDM  and  CDMA  handles  such  channels  [28]: 

“  a  new  waveform  for  UWB  systems  [is]  obtained  by  the  combination  of  Orthogonal  Fre¬ 
quency  Division  Multiplex  (OFDM)  and  Code  Division  Multiple  Access  (CDMA).  The  pro¬ 
posed  system,  called  Spread  Spectrum  -  Multi-Carrier  -  Multiple  Access  (SS-MC-MA)  turns 
out  to  be  a  judicious  solution  to  combat  frequency  selectivity  and  narrowband  interferes,  and 
to  manage  the  coexistence  of  several  users  and  piconets.” 

This  waveform  has  been  under  intense  development  because  of  the  degrees  of  freedom  [14]  and  substantial 
waveform  designs  for  cognitive  radios  [76], 

Regardless  of  the  system,  the  ensemble  of  narrowband  channels  {G(x:  fc)}  stressed  the  distributed 
wireless  power  control  algorithms  and  showed  the  sensitivity  to  system  delay.  Likewise,  the  ensemble  of 
wideband  channels  {G(x:  /)}  should  also  determine  performance  bounds  to  benchmark  these  wideband 
systems. 
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APPENDIX  A 

PROOF  OF  THE  PARETO  FRONT  THEOREM 


Theorem  4  claims  the  Pareto  front  of  the  power-rate  image  is  upper  boundary  of  parameterized  by  the 
mapping 


ps  H >  max 


{-Re(Pc)  :  Pc  €  ^c;ps  =  lTPc}  =:  ^(pe)- 


The  proof  shows  the  mapping  ps  fi(ps)  is  continuous  and  strictly  increasing.  As  such,  the  graph  of  this 
function  is  its  own  Pareto  front  by  Lemma  2.  In  more  detail,  let  pm ;n  :=  lTpmin  and  denote  the  graph  of 
ps  fl(ps)  as 


0(fi)  :  = 


PE 

^(Ps) 


:  PE  A  pm;n 


Let  Q  C  M2  denote  the  quadrant  pointing  to  “less  power  and  more  rate”. 


Q:= 


-PE 

rs 


:  Ps,rs  >  0 


Because  this  graph  f/(Q)  is  its  own  Pareto  front,  the  Contact  Theorem  (Theorem  1)  requires  any  contact  of 
this  quadrant  Q  located  at  a  point  on  the  graph  to  be  precisely  that  point  of  contact: 


Ps 

.^(Ps) 


PS 

.^(Ps) 


+q  n  gw. 


Because  this  upper  boundary  is  parameterized  by  continuous  and  strictly  increasing  function  pe  fl(p s), 

contact  with  the  power-rate  image  (Definition  9) 


m(pc) 


lTPc 

.-Re(Pc). 


:  Pc  €  Pc 


must  still  be  the  single  element: 


PE 

D(ps) 


PE 

.^(ps) 


+q  n 


The  Contact  Theorem  (Theorem  1)  applies  again  to  yield  that  point  on  the  graph  Q  (Q)  must  be  a  Pareto 
point  with  respect  to  the  power-rate  image  ffi(T’c).  Consequently,  the  graph  G(Cl)  is  the  Pareto  front  of 
the  power-rate  image  91  (Vc)-  The  remainder  of  this  appendix  is  devoted  to  showing  O  is  continuous  and 
increasing. 


Lemma  5  Let  X  and  Y  be  compact  topological  spaces.  Let  h  :  X  x  Y  R  be  continuous.  Define 

fl(y)  :=  max{/i(x,y)  :  x  G  X}. 

Then  Q  :  Y  -tlij  well-defined  and  continuous. 


Proof:  Because  h{ o,  y)  is  continuous  on  the  compact  set  X ,  the  maximum  is  well  defined.  Suppose  {yn} 
is  a  sequence  in  Y  converging  y  £  Y.  For  each  n  =  1,  2,  ...  define  xn  £  X  as  a  maximizer:  h(xn,  yn)  = 
D(yn).  By  compactness  of  X,  {x„}  contains  a  convergent  subsequence.  By  relabeling,  let  xn  — >  x  £  X . 
The  corresponding  subsequence  of  {yn},  also  relabeled,  still  converges  as  yn  — >  y  £  Y.  Because  h  is 
continuous  on  X  x  Y, 


lim  h(xn,  y 

n— >oo 


)  =  h(x,y). 


A-l 


Consequently,  the  sequence  (Q(yn)}  is  a  convergent  sequence  with  limit 

lirn  fi(yn)  =  lim  h(xn,  y  )  =  h(x,y). 

n— >00  n— >oo 


By  definition, 

7 

Select  any  x*  G  X  that  is  a  maximizer 
By  definition, 

Taking  the  limit  on  both  side  produces 


lim  fi(yn)  =  h(x,  y)  <  ft(y). 

I — >-oo 

h(x*,y)  =  n(y). 
n(yJ  >  Mx*>yJ- 


lim  0(yn)  >  lim  /i(x*,yn)  =  /i(x*,y) 


^(y). 


Combine  both  inequalities  to  get 
or  that  f2(yn)  converges  to  f2(y).  Ill 


n{y)  <  lim  Q(yn)  <  fi(y) 


Lemma  6  Let pm ;n  :=  lTpmin.  Define  17  :  [pm;n,  oo)  — >  M+  as 

ps  i— >  max  {.Rs(Pc)  :  Pc  G  'PcSPs  =  lTPc}  =:  ^(Ps)- 

Then  Q,  is  continuous. 


Proof:  The  feasible  power  vectors  pc  G  Vc  are  parameterized  as 

PC  =  Pmin  +  (I  ~  i?)'1P A-!  PK  €  K+- 

Let  Uk  C  consisting  of  all  power  vectors  u k  with  unit  network  power: 

1  =  1  t(I-F)-1uK- 

Observe  that  Uk  is  closed  and  bounded  because  (/  —  F)  is  invertible.  Equivalently,  Uk  is  compact.  Let 

Pmax  5“  Pmin-  Define  (j)  [pmin>Pmax]  X  Uk  t  F C  as 

0(ps,ujc)  :=  pmin  +  (ps  -  Pmin) (I  ~  f)_W- 

Because  i?s  °  0  is  continuous  on  the  compact  domain  [pmin,Pmax]  x  Uj<.  Lemma  5  determines  that 

ps  — >•  max  {  /is  ( (p (ps ,  u  c )  :  ux  G  ZCfr} 

must  be  continuous.  Observe  that  (j)  parameterizes  each  ps-slice  of  the  feasible  power  vectors: 

1T</>(PS,  UK)  =  lTpmin  +  (ps  -  Pmin)lT(-f  -  F)-1Vk 
=  Pmin  T  (PE  —  Pmin) 

=  PS 

for  Pmin  <  Ps  <  Pmax-  Consequently, 

O(ps)  =  max  |i?s(Pc)  :  Pc  e  'PciPs  =  lTPc} 

=  max{i?s(0(ps,UA')  :  ur'  G  Z7r'} 

or  that  0[pmjn,  Pma.x]  — >  M+  must  be  continuous.  Because  pmax  is  arbitrary,  0[pmjn,  oo)  — >•  M+  must  also 
be  continuous.  Ill 
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Lemma  7  Let  pm\n  :=  lTpmin.  Define  fl  :  [pmin,  oo)  —>  M+  as 

Pt,  h-  max  |-Rs(Pc)  :  Pc  G  Vc-,Ps  =  lTPc}  =:  ^(pe)- 
Then  O  is  strictly  increasing. 

Proof:  Set  /7v:  >  lTpmin.  Let  pmax  be  any  maximizing  vector  of 

max  |f?s(PcO  :  Pc  e  Vc;ps  =  lTPc}  ■ 

Compactness  of  the  the  intersection  of  the  feasible  power  vectors  Vc  with  the  linear  constraint  coupled 
with  the  continuity  of  R.y  guarantees  existence  of  at  least  one  pmax.  Set  fix  ■=  7iw(pmax)  an^ 

F  =  G^rvG/;  rx  =  diag(7x). 

Because  fixfillx-  the  spectral  radius  of  F  cannot  exceed  the  spectral  radius  of  F.  By  Theorem  3,  the  feasi¬ 
ble  power  vectors  p(  6  Tlfi  exceeding  this  new  SINR  constraint  fixfiTAtiPr)  belong  to  the  cone 

P C  G  Vc  ■■=  PX  +  (I  -  F)-1^ 


with  vertex 

Px  =  (I  ~  Fy'G^TxPN. 

By  Theorem  3,  this  vertex  is  the  unique  solution  to  the  equality  fix  =  7m  (Pc)  f°r  p q  €  Vc-  Thus, 
Pmax  =  Px  and  y.uiPx)  =  lx-  That  is,  7m (Px)<7m(Pc)  with  equality  if  and  only  if  px  =  pr. 
Consequently,  -Re(Px)  A  -^s(Pc)  with  equality  if  and  only  if  px  =  Pc-  Because  p  v  is  the  vertex  of  Vc, 
lTpx  <  lTPc>  with  equality  if  and  only  if  px  =  Pc-  Select  any  A  pc  >  0.  The  preceding  inequalities 
force 

-Rs(Pa')  <  max  {-Rs(Pc)  :  Pc  e  ^c,Ps  +  A pc  =  lTpc}  • 

By  construction,  Vc  C  Vc  so  that 

-Rs(Px)  <  max  j-Rs(Pc)  :  Pc  G  Vc,Pz  +  A pc  =  lTpc}  • 

Consequently,  the  mapping 

Py  max  {-Rs(Pc)  :  Pc  e  VC]Pz  =  lTPc} 

is  a  strictly  increasing  function.  Ill 
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